大 连理 工大 学 本 科 叶 业 论 文 


铅 忆 弹 高 速 冲击 钢 靶 仿真 研究 


Numerical Simulation of Lead Core Projectile Impacting Steel 


Target at High Velocity 
学 Bú (ЖО: 机 械 工 程 学 院 
专 МИ: 机 械 设计 制造 及 其 自动 化 创新 实验 班 
学 生 WE E EAR GH 
ae =: 201962002 
指导 教师 : TEC HR 
评阅 教师 : 35 9] fe 


完成 日 期 : 2023 年 7 月 10 H 


大 连理 工大 学 


Dalian University of Technology 


铅 蕊 弹 高 速 冲击 钢 轩 仿真 研究 


原创 性 声明 


本 人 郑重 声明 : 本 人 所 呈 交 的 毕业 论文 (设计) ， 是 在 指导 老师 
的 指导 下 独立 进行 研究 所 取得 的 成 果 。 毕 业 论文 《设计 ) 中 凡 引 用 他 
人 已 经 发 表 或 未 发 表 的 成 果 、 数 据 、 观 点 等 ， 均 已 明确 注 明 出 处 。 除 
文中 己 经 注 明 引用 的 内 容 外 ， 不 包含 任何 其 他 个 人 或 集体 已 经 发 表 或 
撰写 过 的 科研 成 果 。 对 本 文 的 研究 成 果 做 出 重要 页 献 的 个 人 和 集体， 
均 已 在 文中 以 明确 方式 标明 。 


本 声明 的 法 律 贡 任 由 本 人 承担 。 


作者 签名 Fath Н Hj. 20234E 7 H 10 H 


A LCS E A АҒА, 


关于 使 用 授权 的 声明 


本 人 在 指导 老师 指导 下 所 完成 的 毕业 论文 (设计) 及 相关 的 资料 
(包括 图 纸 、 试 验 记 录 、 原 始 数 据 、 实 物 照片 、 图 片 、 录 音 带 、 设 计 
手稿 等 ) ， 知 识 产 权 归 属 大 连理 工大 学 。 本 人 完全 了 解 大 连理 工大 学 
有 关 保 存 、 使 用 毕业 论文 (设计) 的 规定 ， 本 人 授权 大 连理 工大 学 可 
以 将 本 毕业 论文 《设计 ) 的 全 部 或 部 分 内 容 编 入 有 关 数 据 库 进行 检索 ， 
可 以 采用 任何 复制 手段 保存 和 汇编 本 毕业 论文 (设计) 。 如 宁 发 表 相 
天成 果 ， 一 定 征 得 指导 教师 同意 ， 且 第 一 筑 名 单位 为 大 连理 工大 学 。 
本 人 离 校 后 使 用 毕业 毕业 论文 (设计) 或 与 该 论文 直接 相关 的 学 术 论 
文 或 成 果 时 ， 第 一 署名 单位 仍然 为 大 连理 工大 学 。 


论文 作者 签名 : FRX H Hj: 20234E 7 H 10 H 
әже. B BA As Н 期 20234 7 H 10 H 


EE ТА A ҒА, 


摘 =< 


高 速 侵 彻 数值 模拟 研究 通常 采用 拉 格 朗 日 有 限 元 法 。 < 
Вр Е 2. 适当 的 
2... 72 
流体 动力 学 (SPH) 方法 模拟 5.56 mm M193 步枪 弹 侵 彻 钢板 过 程 ， 研 究 单 元 类 型 与 仿 
真 算法 对 模拟 精度 的 影响 ,分析 铬 芯 弹 穿 甲 过 程 和 原理 ,并 与 实验 结果 比较 。 结 果 表 明 ， 
SPH 法 在 模拟 铅 芯 弹 高 速 冲 击 问 题 中 具有 最 高 的 仿真 精度 、 稳 定性 和 效率 ， 可 准确 复 现 
实验 中 靶 板 的 变形 、 侵 蚀 、 剖 塞 过 程 。 高 速 铅 心 弹 对 钢 划 的 破坏 主要 分 为 半 流 体 侵 彻 、 
背 表面 损伤 、 剪 切 冲 塞 三 部 分 : 400 HBW 及 以 下 的 中 低 硬 度 装甲 钢 以 前 两 者 为 主 ，400 
HBW 部 由 于 拉 伸 塑性 较 低 而 易 发 生 背 面 月 落 ; 500 HBW 及 以 上 的 高 硬度 装甲 钢 以 后 者 
为 主 ， 界 面 击溃 发 生 于 500-600 HBW 间 。 提 出 基于 Rosenberg 半 流 体 侵 彻 理论 和 Recht 
I Eee \ 弹 穿 甲 模型 ， 受 限于 模型 本 身 的 不 准确 性 和 弹 体 几何 的 复杂 性 ， 仅 
能 定性 描述 穿 甲 原理 ， 以 及 初步 预测 铅 芯 弹 对 靶 板 的 破坏 模式 。 


关键 词 : 光滑 粒子 流体 动力 学 ; 有 限 元 仿真 ; 高 速 冲击 ; 半 流 体 侵 彻 
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Numerical Simulation of Lead Core Projectile Impacting Steel Target 
at High Velocity 


Abstract 


Numerical simulation of high-velocity penetration is usually carried out by Lagrangian 
finite element method. The process of lead-core full metal jacket bullet impacting steel target 
belongs to a typical large deformation and high strain rate problem, and the traditional 
Lagrangian method has poor simulation effect on it. In order to find a suitable simulation 
method, this paper adopts Lagrangian finite element method (FEM), element-free Galerkin 
(EFG) and smoothed particle hydrodynamics (SPH) to simulate the process of 5.56 mm M193 
penetrating steel target, studies the influence of element form and simulation algorithm on the 
accuracy, analyses the penetration process and principle of lead-core bullet, and compares 
with the experimental results. Results show that SPH has the best accuracy, stability and 
efficiency in simulating the high-velocity impact problem of lead-core bullet, and can 
accurately reproduce the deformation, erosion and plugging process of the target plate in 
experiment. The damage of steel target by high-velocity lead-core bullet can be divided into 
three parts: semi-fluid penetration, back-surface chipping and shear plugging. The former two 
are dominant for medium and low hardness armour steel of 400 HBW and below, and the 
back-surface chipping is easy to occur in the target of 400 HBW due to its low tensile 
plasticity; while the latter is dominant for high hardness armour steel. A penetration model of 
lead-core bullet based on Rosenberg semi-fluid penetration theory and Recht shear plugging 
theory is proposed. Limited by the inaccuracy of the model itself and the complexity of the 
projectile geometry, it can only describe the penetration principle qualitatively and predict the 
damage pattern of lead-core bullets on the target plate preliminarily. 


Key Words: Smoothed Particle Hydrodynamics; Finite Element Method; High-velocity 


Impact; Semi-fluid penetration 
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铅 芯 一 铀 被 弹 〈Lead-core Full Metal Jacket) 是 最 典型 的 一 种 步枪 弹 结 构 ， 其 起 源 可 
追溯 至 19 世纪 末期 的 瑞士 的 施 密 特 一 鲁 宾 步枪 ， 旨 在 提高 枪弹 整体 强度 以 减少 在 内 弹 
道 过 程 中 的 弹 体 变形 和 枪 管 磨损 。 随 着 单 兵 及 轻型 车 辆 防护 技术 的 发 展 ，FMJ FPR 
SIN A ACH, UE AM eR, т, HP AAA RIP! Spuren 
度 高 站 的 优点 ， 时 至 今日 ，FMJ 结构 仍 被 广泛 应 用 于 下 人 至 小 口径 步枪 、 上 至 高 射 机 枪 的 
各 类 步 /机 枪弹 中 。 除 此 以 外 ， 铅 芯 弹 穿 甲 性 能 较 钢 一 铅 芯 组 合 的 半 穿 甲 弹 稳定 内， 因此 
也 被 大 量 用 于 评估 轻型 装甲 的 抗 弹 能 力 口 。 

要 对 目标 造成 有 效 破坏 ， 枪 弹 必 须 具 有 足够 的 威力 以 击 穿 相应 的 防护 装备 。 铅 蕊 弹 
的 主要 杀伤 对 象 为 装备 有 防弹 衣 等 的 有 生 目 标 , 其 防护 设备 材料 通常 为 钢 、 纤维 和 陶 次 。 
早期 单 兵 防护 设备 通常 由 高 强度 钢板 冲压 而 成 ， 厚 度 一 般 在 2 mm 以 下 ， 仪 能 有 效 阻挡 
手枪 弹 。 随 着 成 分 设计 、 冶 炼 工艺 和 热处理 工艺 的 不 断 发 展 ， 高 性 能 装甲 钢 不 断 涌现 。 
从 中 低 强 度 向 高 硬度 、 超 高 硬度 、 高 强度 和 高 韧性 发 展 铅 。 当 代 装 甲 钢 依 和 硬度 分 为 均 质 
轧 制 装甲 (Rolled Homogeneous Armour) 、 高 硬度 装甲 钢 CHigh-Hardness Armour) 和 
超 高 硬度 装甲 钢 C Ultra-High-Hardness Armour) 三 类 。RHA 硬度 在 400 HBW 以 下 ,在 
保证 综合 抗 弹 能 力 的 同时 兼 具 良好 的 焊接 、 冷 弯 和 切割 性 能 ， 主 要 用 于 制造 几何 形状 复 
杂 且 厚度 较 大 的 坦克 装甲 四，HHA 硬度 约 为 500 HBW， 兼 顾 硬度 和 韧性 ， 广 泛 用 于 制 
造 各 类 单 兵 防护 设备 〈 如 防弹 衣 ) 由 和 坦克 、 装 甲 车 等 的 复合 装甲 ，UHHA 硬度 在 600 
HBW 以 上 ， 多 用 于 轻型 装甲 车 的 披挂 装甲 。 为 综合 体现 铅 心 弹 对 各 类 目标 的 穿 透 能 
ASC AT HE A HERE FE 300~600 HBW |8], 2 zi SLABS FA a AY) E SE FE Y FR] o 

目前 枪弹 穿 甲 研究 以 高 硬度 钢 芯 的 穿甲弹 为 主 ， 对 铅 芯 弹 的 探讨 相对 较 少 ， 且 多 局 
限于 实验 。 Kiis 等 外 通过 实验 与 数值 模拟 研究 了 7.62 mm 钢 蕊 穿甲弹 对 高 强度 钢 靶 的 侵 
TIGER. AM SPH 模拟 的 靶 板 对 弹 体 的 阻力 较 FEM 法 偏 低 ， 但 能 更 好 地 表现 崩落 和 碎 
片 飞溅 过 程 。 然 而 SPH 存在 多 种 算法 ， 以 LS-DYNA 为 例 ， 其 可 选 近似 方法 总 计 多 达 
16 种 ， 而 原文 仅 使 用 默认 方法 ， 未 做 进一步 探究 。 李 勇 1" 通 过 数值 模拟 研究 14.5 mm 
钨 合 金 罕 甲 弹 对 钢 部 的 侵 彻 过程 ， 发 现 当 使 用 对 称 边 界 时 ，FEM、SPH、SPH-FEM ЖҚ 
合 三 种 方法 计算 结果 相差 很 小 ， 均 能 准确 复 现 延 性 穿 甲 产 生 的 花 浙 形 翻 唇 ;SPH 的 计算 
耗 时 高 于 FEM。Bervik 等 U1 研究 了 7.62 mm 铅 芯 弹 和 钢 芯 弹 对 五 种 高 强度 钢 的 侵 彻 过 
程 ， 发 现 单 屋 和 无 间 辽 的 双 层 靶 具 有 相同 的 抗 弹 性 能 ， 而 有 间 隐 的 双 层 靶 劣 于 前 两 者 ; 

维 轴 对 称 的 FEM 和 多 物质 ALE 都 不 能 可 靠 地 模拟 铅 芯 弹 的 侵 彻 过 程 。Magier 51121 
使 用 单 物质 ALE-FEM 耦合 法 模拟 7.92 mm DS 反 坦 克 步 枪弹 对 早期 装甲 钢 的 侵 彻 过 程 ， 
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与 历史 数据 对 照 具 有 良好 的 精度 ;并 发 现 高 速 下 铅 世 与 靶 板 会 相互 侵蚀 ， 铬 芯 与 被 甲 的 
碎片 沿 孔 壁 反 向 流动 ， 过 程 类 似 于 高 速 长 杆 弹 的 半 流 体 侵 彻 。Coghe 等 5 通过 实验 获得 
T 5.56 mm M193 铅 芯 步 检 弹 在 不 同 速度 下 对 6.5 mm Secure 500 的 毁伤 形 貌 ,发现 铅 蕊 
弹 在 高 速 下 对 高 硬度 装甲 钢 的 破坏 模式 以 剪 切 冲 塞 为 主 ， 表 面 仅 发 生 少 量 侵蚀 ， 且 穿 甲 
能 力 较 前 钢 芯 、 后 铅 芯 的 55109 更 优秀 ，Sherrill 等 0 通过 实验 发 现 较 低 速度 下 两 种 枪 
弹 表 现 与 高 速 下 相反 ， 且 铅 芯 弹 对 陶瓷 一 纤维 复合 装甲 穿 甲 效 能 显著 劣 于 半 钢 芯 弹 。 
Binar 等 03 通 过 实验 获得 了 M193 对 6 mm Armox 600 超 高 硬度 装甲 钢 的 弹道 极限 速度 。 
Siriphala 等 09 通 过 实验 获得 了 M193 对 7 mm Armox 600 的 毁伤 形 貌 ,并 发 现 FEM 无 法 
准确 模拟 侵 彻 过 程 。 可 见 ， 对 铅 心 弹 高 速 穿 甲 的 研究 仍 处 于 初步 阶段 ， 不 同 速度 、 不 同 
误 板 强度 仿真 与 实验 所 得 结论 有 很 大 差别 ， 暂 不 存在 统一 、 广 泛 适 用 的 理论 模型 。 

在 基于 拉 格 朗 日 坐标 的 有 限 元 方法 (finite element method) 中 ， 计 算 网 格 的 每 个 节 
点 在 其 运动 过 程 中 跟随 指定 的 粒子 ， 因 此 外 表面 和 接触 面 很 容易 被 追踪 。 由 于 标准 有 限 
元 方法 的 接触 算法 能 可 靠 定 义 ， 且 计算 较为 快速 ， 因 此 是 研究 侵 彻 问题 的 主要 数值 模拟 
工具 。 尽 管 具 有 以 上 优点 ， 拉 格 朗 日 有 限 元 法 对 大 变形 、 大 失效 工 况 一 一 如 本 文 研究 的 
高 速 铅 芯 弹 冲击 问题 一 一 处 理 不 佳 Q0023。 若 网 格 发 生 严重 扭曲 、 变 形 ， 则 可 能 导致 剪 切 
BUEU, 或 由 于 单元 自身 节点 与 面相 互 穿 透 而 引发 负 体积 错误 ， 进 而 使 单元 在 发 生物 理 
失效 前 先 数值 失效 ， 以 致 抗 弹性 能 异常 偏 低 。 此 外 ， 由 于 网 格 法 必须 删除 失效 单元 ， 整 
个 计算 过 程 中 质量 和 动量 不 守恒 ， 故 不 适用 于 对 弹 诅 惯性 敏感 的 高 速 冲击 问题 08， НІҢ 
于 接触 界面 的 不 断 变化 而 引发 非 真实 的 接触 力 振荡 021， 进 而 可 能 导致 弹 靶 过 度 侵蚀 。 此 
外 ， 对 于 高 速 侵 彻 问题 ，FEM 还 存在 部 件 间接 触 穿 透 Pa 和 接触 力 偏 低 避 的 问题 。 无 网 
格 伽 辽 金 方法 (Element-Free Gelarkin) 部 分 摆脱 了 对 网 格 的 依赖 。 由 于 节点 间 不 存在 连 
接 ， 避 免 了 单元 扭曲 问题 ， 较 FEM 显著 提升 了 计算 稳定 性 。 然 而 尽管 EFG 被 归 类 为 无 
网 格 方法 , 但 依赖 背景 网 格 以 进行 数值 积分 , 因此 仍 存在 由 网 格 /单元 删除 引发 的 质量 不 
守恒 等 问题 。 与 其 它 无 网 格 方法 类 似 ，EFG 需要 对 边界 条 件 进行 特殊 处 理 24， 其 中 由 于 
与 标准 拉 格 朗 日 有 限 元 法 的 相似 性 ，FEM-EFG 耦合 被 认为 是 一 种 有 效 的 方法 。 

除 拉 格 朗 日 坐标 外 ， 欧 拉 坐 标 类 方法 由 于 对 大 变形 的 优秀 适应 性 也 在 数值 仿真 中 受 
到 广泛 应 用 。 目 前 在 高 速 侵 彻 问题 中 主要 应 用 的 是 任意 拉 格 朗 日 一 欧 拉 方法 (Arbitrary 
Lagrange-Euler) 。 但 ALE 法 的 材料 边界 不 清晰 ， 对 强度 模拟 也 不 准确 中 3， 因此 仅 适 用 
于 强度 较 小 且 大 变形 的 枪弹 ， 而 不 可 用 于 模拟 高 强度 、 小 变形 的 鞭 板 ， 故 必须 使 用 流 固 
#44 (Fluid-Solid Interaction) 以 模拟 两 类 材料 的 交互 。 然 而 ， 所 有 ALE 方法 一 一 包括 
单 物 质 ALE、 多 物质 ALE 和 结构 化 ALE 一 一 在 与 固体 单元 耦合 时 均 存 在 流体 泄露 问题 ， 
必须 依 经 验 设置 参数 才 可 确保 稳定 多 ， 无 法 被 用 于 工 况 复杂 多 变 的 大 规模 分 析 。 此 外 ， 
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由 于 ALE 是 Euler 方法 的 一 种 变 体 ， 其 单元 和 节点 在 整个 模拟 过 程 中 固定 、 材 料 在 共 节 
点 的 单元 间 流 动 ， 而 非 如 拉 格 朗 日 类 方法 一 般 网 格 和 材料 协同 运动 ， 无 法 直接 由 节点 准 
确 描述 部 件 的 形状 。 
光滑 粒子 流体 动力 学 法 (Smooth Particle Hydrodynamics) 是 一 种 拉 格 明日 无 网 格 方 
法 。 通 过 一 组 插值 粒子 离散 系统 体积 ， 每 个 粒子 携带 着 计算 场 和 内 在 数据 的 离散 值 ， 并 
在 控制 方程 的 规律 下 移动 。SPH 作为 一 种 完全 的 无 网 格 方法 ， 能 从 根本 上 避免 网 格 扭曲 
带 来 的 一 系列 问题 ， 同 时 ， 由 于 采用 Lagrange 坐标 ， 也 不 存在 Euler 类 方法 难 定 位 的 缺 
点 。 因 此 ， 从 以 上 来 看 ，SPH 是 一 种 适用 于 模拟 铅 芯 弹 冲 击 钢 靶 的 方法 。 然 而 ， 更 广泛 
的 研究 表明 ，SPH 同样 存在 根本 性 缺陷 ,在 高 速 冲 击 中 主要 表现 为 应 力 不 稳 定 和 粒子 异 
常 穿 透 。 应 力 不 稳 定性 的 根本 原因 在 于 应 力 状 态 与 核 函数 的 不 匹配 : 负 压 状态 下 粒子 间 
产生 随 着 粒子 间距 的 减 小 而 增 大 的 引力 ， 导 致 拉 伸 不 稳定 性 ; 正 压 状态 下 粒子 间 产 生 随 
着 粒子 间距 的 减 小 而 先 增 大 后 减 小 的 斥 力 ， 导 致 压缩 不 稳定 性 。 杨 秀峰 等 上 5 改进 了 传统 
SPH 方法 的 核 函 数 和 离散 格式 , 从 原理 上 消除 了 此 问题 , 并 具有 理论 上 的 普 适 性 。 然 而 ， 
其 改进 方案 目前 尚未 经 过 工程 验证 ， 暂 无 法 确定 是 否 存在 其 它 方面 的 缺陷 。Xnu 等 eq 在 
LS-DYNA 的 基础 上 开发 了 自 适 应 、 重 归 一 化 和 Lagrange 近似 方法 ， 在 配合 四 阶 样 条 核 
函数 时 ， 三 者 均 能 较 好 地 抑制 了 拉 伸 不 稳定 的 问题 。Lagrange 近似 法 仅 可 被 用 于 求解 小 
变形 问题 ， 自 适应 核 函数 需 搭 配 使 用 椭圆 支撑 域 ， 其 计算 稳定 性 不 及 传统 球形 支撑 域 ， 
在 发 生 大 量 失 效 时 此 缺陷 尤为 明显 ; 重 归 一 化 近似 广泛 适用 于 绝 大 多 数 固体 问题 ,但 对 
拉 伸 不 稳定 的 抑制 效果 不 及 前 两 种 方法 。Liu 等 R71 使 用 重 归 一 化 近似 方法 模拟 了 有 限 厚 
和 半 无 限 混凝土 对 受 尖 头 弹 侵 彻 的 过 程 , 准确 复 现 了 实验 中 靶 板 正面 变形 与 靶 后 碎片 云 
的 形 貌 。Bojanowski 等 PH 研究 长 杆 弹 侵 彻 土壤 的 过 程 ， 发 现 SPH 和 多 物质 ALE 均 能 获 
得 准确 的 模仿 真 结果 ， 其 中 SPH 能 与 FEM 良好 耦合 和 接触 ， 并 且 对 网 格 密度 和 时 间 步 
长 不 敏感 。 然 而 ， 以 上 对 SPH 的 应 用 研究 多 局 限于 刚体 侵 彻 ， 弹 体 使 用 无 变形 和 失效 
BJ FEM 网 格 ， 其 结论 未 必 完 全 适用 于 接触 面 两 侧 均 存 在 大 量 失效 的 高 速 半 流 体 侵 彻 。 
综 上 所 述 ， 铬 芯 弹 高 速 冲击 /问题 在 理论 和 仿真 上 均 不 成 熟 。 本 文 利用 LS-DYNA s 
现 FEM、EFG、SPH 三 种 方法 模拟 5.56X45 mm M193 步枪 弹 冲 击 钢 贰 的 过 程 ， 探 究 单 
元 算法 与 参数 设置 对 仿真 精度 的 影响 , 基于 以 上 所 得 模型 研究 铅 蕊 弹 侵 彻 不 同 强 度 钢 误 
板 的 过 程 ， 参 考 半 流体 侵 彻 理论 和 仿真 结果 提出 铅 芯 弹 穿 甲 公式 ， 为 揭示 铅 忌 弹 穿 甲 机 
里 和 针对 性 防护 设计 优化 提供 一 定 借 鉴 思 路 。 


ТЕ 
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1 数值 模型 建立 


1.1 有 限 元 模型 


由 于 SPH 方法 对 边界 处 理 不 佳 ， 且 计算 耗 时 远 高 于 传统 网 格 法 PCJ， 因 此 本 文中 所 有 
SPH 模型 实际 均 为 SPH-FEM 耦合 模型 , 仅 大 变形 .大 失效 的 弹 体 及 靶 板 中 心 区 采用 SPH 
粒子 单元 ; 变形 相对 较 小 ， 且 需 设 固定 边界 约束 的 靶 板 外 周 部 分 采用 六 面体 FEM 单元 。 
EFG 及 FEM 弹 体 采用 四 面体 以 准确 描述 其 较 复 杂 的 几何 外 形 和 避免 负 体积 错误 ， 部 板 
整体 采用 六 面体 单元 以 提高 计算 精度 。 所 有 模型 靶 板 中 心 区 及 弹 体 原始 网 格 边 长 均 为 
0.25 mm， 靶 板 外 周 网 格 由 内 向 外 逐渐 稀 疏 ， 最 大 边 长 不 超过 5 mm. SPH 部 件 间 使 用 粒 
子 近似 方法 自动 识别 接触 CCONT=0) ， 无 需 额 外 定义 ; SPH 靶 板 核心 区 SPH 粒子 与 外 
周 区 六 面体 单元 设置 固 连 点 面 接触 (CONTACT TIED NODES TO SURFACE) ， 其 中 中 
心 区 为 从 节点 ， 外 周 区 为 主 面 。FEM、EFG 将 弹 、 部 全 体 部 件 设 为 一 部 件 集 (SET PART 
LIST) ， 定 义 侵 蚀 单 面 接触 (CONTACT ERODING SINGLE SURFACE) ， 此 接触 方法 
会 自动 识别 所 有 接触 面 ， 避 免 碎 片 与 母体 相互 穿 透 。 此 外 ， 设 置 内 接触 CONTACT 
INTERIOR) 以 抑制 网 格 内 穿 透 。 接 触 搜索 深度 DEPTH 设置 为 5， 面 一 面 与 边 一 边 穿 透 
在 计算 时 均 被 检查 。 其 余 参 数 保持 默认 。 

SPH 粒子 近似 方法 采用 重 整 化 近似 (FORM=1) ， 此 方法 适用 于 大 多 数 固体 问题 ; 
核 函 数 设 为 四 阶 样 条 曲线 (ӨРНКЕКМ-І), ЖЕЛ CCSLHO 初 设 为 1.3， 此 组 合 能 
较 好 地 抑制 拉 伸 不 稳定 问题 ,避免 SPH 靶 板 核心 区 在 与 FEM 外 周 区 固 连 的 边界 发 生前 
切 数值 失效 。FEM 六 面体 单元 均 选 择 八 节 点 常 应 力 单元 (ELFORM=1) ， 此 单元 较 全 
积分 单元 (ELFORM=2、-1、-2) 能 更 好 地 避免 因 负 体积 而 导致 的 数值 失效 ， 但 因 是 缩 
减 积分 而 存在 沙漏 问题 ， 故 在 沙漏 控制 (CONTROL HOURGLASS) 中 设置 体积 精确 的 
Flanagan-Belytschko 刚性 沙漏 控制 СІНО-5). EFG 单元 均 选 择 标准 无 网 格 实体 函数 

(ELFORM=41) ， 核 函数 的 归 一 化 扩张 参数 DX、DY、DZ) 均 设 为 1.05， 边 界 处 理 
方法 在 两 组 模型 中 分 别 选 择 FEM-EFG 耦合 法 СЕВТ-3) 与 快速 转换 法 (IEBT=4) , 
其 中 前 者 在 边界 处 理 上 较 可 靠 ， 后 者 具有 最 优 的 计算 效率 和 稳定 性 ， 分 别 记 为 EFG3、 
EFG4; 领域 积分 方法 选择 改进 的 高 斯 积分 CIDIM=3) 。 


(a) SPH 弹 体 (b) FEM/EFG 弹 体 
图 1.1 弹 体 模型 
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(а) SPH-FEM 3E (b) Z FEM/FEM-EFG 3E 
图 1.2 靶 板 模型 


1.2 材料 模型 

弹头 及 靶 板 材料 均 Johnson-Cook 本 构 模 型 和 Gruneisen 状态 方程 。Johnson-Cook 模 
型 广泛 应 用 于 机 械 加 工 、 爆 炸 与 高 速 冲 击 领 域 ,， 综合 考虑 了 材料 的 应 变 强 化 效应 和 应 变 
率 与 温度 对 强度 与 塑性 的 影响 。 其 中 应 力 一 应 变 关系 的 展开 形式 为 : 


ж 


E T-T 
o, =(A+Be" )(1+ Cln —)[1- (一 一 一 )” 1.1 
«7 (cC Be Ee Сз 291-0070 (1.1) 


AP: А. B. n. C. m 均 为 材料 参数 ， 由 实验 确定 。4 为 准 静态 下 屈服 应 力 ; В 为 应 变 
硬化 常数 ，n 为 应 变 硬化 指数 ; C 为 应 变 率 硬化 系数 ; m 为 温度 软化 指数 。cer、se 分 别 
为 等 效应 力 和 应 变 ; e^. eV 分 别 为 实时 和 参考 应 变 率 ; T. Т, 公 分 别 为 实时 、 参 考 和 熔 
Mi BE, Т.Н 22 ‘C， 各 种 钢 熔 点 统一 取 1537.9 С. 

ANE e, 由 式 (1.2) 确 定 : 
E T-T 


T Ет (1.2) 


т ғ 


e, = (D, + D, expD,o`)(1+ D, In 


式 中 : Di Ds 为 材料 参数 ， 由 实验 确定 。 с*= 23 为 应 力 三 轴 度 ， 其 中 ar 为 静水 压力 。 


当 损 伤 变量 D 达到 1 时 单元 发 生 失效 ， 其 计算 方法 为 ; 


E 了 Ағ, 
У. 


(1.3) 


E 


ЖОҒА E Xd 
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表 1.1 Johnson-Cook 强度 模型 主要 参数 
材料 A/MPa B/MPa n C m 
CuZn30P5! 111.7 504.7 0.42 0.009 1.68 
4g Pl 1 63 0.163 0.41 1 
43401301511 719 456 0.093 0.008 1.03 
Mars 19071 1193 500 0.676 0.00435 1.17 
Armox 50013317 1470 702 0.199 0.00549 0.811 
Armox 6008 1580 958 0.175 0.00877 0.712 
表 1.2 Johnson-Cook 失效 模型 主要 参数 
材料 Di D» Ds D4 Ds 
CuZn30D8] 0 2.65 -0.62 0.28 0 
4g P9] 0 0 0 0 0 
43408081 | -0.8 2.1 -0.5 0.002 0.61 
Mars 190821 0.21 7.2 -5.44 0 0 
Armox 5003] 0.068 5.382 -2.554 0 0 
Armox 60081 -0.4 1.5 -0.5 0 0 


1.3 仿真 可 靠 性 验证 


实验 表明 , 6 mm Armox 600 fon] (ul 


阻挡 937 m/s 的 M193P1, 弹道 极限 速度 (V50) 


约 为 990 m/sU?l, HF Armox 600 钢板 存在 0.0/+0.6 mm 制造 公差 63， 故 可 认为 当 子 弹 
速度 为 990 m/s, Es 27 6.0 mm 时 能 恰好 击 穿 , 剩余 速度 、 动 能 接近 于 0。 采 用 FEM. 
EFG3、EFG4、SPH 四 种 方法 模拟 以 上 两 种 速度 工 况 ， 残 余 质 量 (m. . WE (РО. 

动能 (Ew) 、 开 坑 深 度 Dy 与 最 大 直径 Dmax WK 1.3 所 示 。 其 中 SPH 法 穿 出 体 由 靶 板 与 


弹 体 的 碎片 组 成 ， 两 部 分 粒子 的 速度 差异 很 大 。 昌 


于 穿 出 体 以 塞 块 为 主 ， 故 速度 由 其 顶 


部 上 随机 4 节点 轴 疝 速度 取 平 均值 近似 求 得 。 因 仅 有 少量 弹 体 穿 出 ， 而 反 向 、 径 向 飞散 


的 弹 体 碎片 也 具有 


定 动 能 ， 知 统计 在 内 会 产生 很 大 误差 。 因 此 残余 动能 以 丢 体 动能 近 


似 估 计 。 考 虑 到 靶 体 的 整体 运动 ，Ew 实际 值 应 当 低 于 近似 结果 。 残 余 质量 以 残余 动能 
与 速度 估算 。SPH 方法 所 得 靶 板 抗 弹 能 力 最 低 ， 击 穿 、 阻 挡 情 况 与 实验 相同 053， 毁 伤 形 
貌 与 实验 结果 最 接近 MI， 邯 板 受 冲 塞 罕 孔 破坏 ， 冲 击 面 仪 受 轻 微 侵蚀 。FEM 法 所 得 击 


穿 、 阻 挡 情况 与 实验 结果 相同 ， 毁 伤 形 貌 与 实验 
锥 形 申 坑 ， 四 坑 底 部 发 生 延 性 穿孔 。 两 种 EFG 法 所 得 靶 板 抗 弹 能 
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ARIZA: 靶 板 受 明 显 侵蚀 ， 形 成 


最 高 ， 击 穿 、 阻 挡 
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情况 与 实验 结果 不 一 致 ， 毁 伤 形 貌 与 实验 结果 相差 较 大 大 : 靶 板 受 明 显 侵蚀 ，EFG3 JÉ 
成 半 柱 形 一 半 锥 形 凹 坑 ，EFG4 НӘР. 


表 1.3 四 种 仿真 方法 比较 
仿真 方法 vms! m/g  v/ms! EW] Рут Dymax/mm 
FEM 97 REF, ФЕН ЛАН 7.1 7.5 
FEM 990 0.04 241 2.51 ee 7.3 
EFG3 97 REF, Wd 4.1 7.0 


SI 


EFG3 990 EF, FRA MIE 6.1 6.8 
EFG4 937 REF, AERA cnt 2.8 7.9 
EFG4 990 REF, MRAM 4.1 7.8 
SPH 97 REF, FU AC 3.1 11.3 
SPH 990 7.0 186 121 He 9.6 


(a) FEM (b) EFG3 (c) EFG4 (d) SPH 
图 1.3 990 m/s M193 4244] 6 mm Armox 600 终 时 刻 弹 靶 Johnson-Cook 损伤 


对 比 四 种 仿真 模型 在 990 m/s 时 靶 板 毁伤 形 钥 与 实验 参考 ， 可 见 以 Armox 600 为 代 
表 的 超 高 硬度 装甲 钢 在 受 铅 忌 弹 冲 击 时 ， 靶 板 损 伤 分 为 挤 压 侵蚀 和 剪 切 冲 塞 两 部 分 ， 实 
验 和 SPH LSI THSE NE. SPH 粒子 在 达到 物理 失效 后 应 力 被 归 零 并 失去 体积 响应 ， 
但 仍 存在 于 积分 域 中 ， 且 保持 与 未 失效 粒子 的 接触 。 失 效 粒 子 在 残余 弹 体 与 训 板 间 形 成 
了 一 层 无 应 力 、 几 平 无 厚度 的 缓冲 层 ， 避 免 了 因 单 元 删除 而 产生 的 非 真实 空 际 ， 单 元 应 
力 不 会 因 接 触 的 突然 断 开 和 恢复 而 振荡 。SPH 粒子 间 的 自动 接触 较 网 格 单元 间 的 接触 更 
平滑 ， 接 触 面 上 压力 更 均匀 ， 减 少 因 非 真实 集中 应 力 导 致 异常 失效 。 以 上 两 点 使 SPH 
鄞 受 挤 压 侵 蚀 更 少 。FEM 及 EFG 法 计算 受 限于 网 格 ， 失 效 单元 在 仿真 过 程 中 被 删除 ， 
因而 不 存在 无 应 力 缓冲 层 ， 邮 板 的 抗 挤 压 侵蚀 能 力 劣 于 实际 情况 。 与 以 压缩 破坏 为 主 的 
长 杆 弹 侵 彻 土壤 1 类 似 , EFG 靶 抗 压 优 于 FEM， 其 中 EFG3 由 于 单元 行为 类 似 于 FEM, 
其 挤 压 破 坏 程 度 也 介 于 FEM 与 EFG4 之 间 。 


E 
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从 图 1.4 阻力 变化 过 程 可 以 看 出 ， 与 前 文 对 靶 板 侵蚀 分 析 类 似 ，SPH 粒子 失效 前 后 
速度 和 动量 不 变 B9， 始 终 保持 对 靶 板 的 神 击 ,因此 铅 芯 和 整个 弹 体 所 受阻 力 均 明 显 高 于 
FEM 和 EFG。 在 冲击 面 直径 与 其 它 方法 相差 不 大 的 前 提 下 ，SPH 部 的 环形 剪 切 面 分 配 
到 更 大 的 剪 应力。 由 于 侵 彻 过 程 中 弹 体 截面 积 和 肢 板 剩余 厚度 均 会 不 断 变 化 ， 无 法 直接 


200 4 


н 
a 
e 


Resistance (kN) 
8 


确定 其 每 个 时 刻 的 具体 值 ， 故 表 14 取 阻 力 峰 值 时 刻 近似 作为 最 大 前 应 力 时 刻 分 析 ， 其 
中 玉 为 最 大 阻力 ， 记 为 实时 接触 面 直径 ，! 为 靶 板 剩余 厚度 ，r 为 实时 型 


sen EFG4 Core 


БЕРЕРІ EFG4 Total 


Ep I Rr RR as Saam TE RR E RR KL I ERU 


Time (ps) 


1.4 990 m/s M193 侵 彻 6 mm Armox 600 过 程 阻 力 


表 1.4 990 m/s БИО Н fovit SE S UTR 


仿真 方法 ЕМ D/mm ¿mm т/МРа 
FEM 55 6.8 5.2 495 


EFG3 71 7.4 3.7 825 
EFG4 85 7.6 4.1 868 


SPH 209 8.6 49 1579 


Не. 


根据 Johnson-Cook 模型 , 在 常温 准 静 态 条 件 下 , Armox 600 的 抗 剪 强度 为 1227 MPa, 


FEM、EFG3、EFG4 所 得 最 大 到 


最 大 到 


论 剪 应 力 均 显 著 低 于 此 值 ， 故 不 发 生 剪 切 破坏 ; SPH 的 


LE 论 均 应 力 高 于 此 值 ， 剪 切面 撕 裂 并 形成 近似 柱 形 塞 块 。 理 论 分 析 与 图 1.3 中 的 念 


真 结果 相 比 可 见 ，EFG3、EFG4 和 SPH 所 得 的 剪 切面 均 不 是 假设 的 规则 圆柱 面 ， 而 是 


受 背 表面 弯曲 产生 的 拉 应 力 与 正 表面 冲击 产生 的 压 应 力 共同 影响 、 直 径 先 增 大 后 减 小 的 
双 锥 组 合 面 ， 其 挤 压 端 直径 略 大 于 弹 坑 直径 。 但 从 图 1.6 (a) 中 的 Johnson-Cook 损伤 和 图 
1.3 (d) 的 最 终 冲 塞 孔 来 看 ， 综 合 应 力 与 应 力 三 轴 度 的 影响 ， 实 际 受 破坏 的 仍 是 与 理论 假 
设 尺 寸 、 形 状 均 相 近 的 圆柱 面 。FEM 法 挤 压 面 侵蚀 严重 ， 呈 近似 半球 形 ， 受 压 区 沿 径 
向 扩张 明显 ， 因 此 剪 切面 仅 有 锥 形 部 分 。 
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Triaxiality Factor (-plvm) 
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图 1.5 990 m/s M193 冲击 6 mm Armox 600 最 大 阻力 时 刻 靶 板 应 力 三 轴 度 


高 速 冲 击 下 郭 板 受 应 变 率 强 化 和 温度 软化 效应 影响 ， 实 际 抗 剪 应 力 随 之 增 减 。 以 与 
实验 结果 最 接近 的 SPH 法 为 例 〈 图 1.6 (Ь)) ， 剪 切面 上 的 Tresca 应 力 约 为 1271 MPa, 
表明 应 变 率 强化 影响 略 大 ， 但 与 准 静态 一 常温 值 相 比 增 幅 很 小 ， 初 步 分 析 时 仍 可 用 后 者 
近似 估算 。 虽 然 理论 最 大 剪 应 力 可 用 于 粗略 判断 是 否 发 生 冲 塞 ， 但 由 于 剪 切 破坏 并 非 仅 
与 强度 有 关 的 瞬间 过 程 ， 因 此 分 析 准 确 性 较 差 。 此 外 ，FEM 和 EFG 虽 在 较 小 变形 情况 
下 计算 效率 较 高 CJ， 但 由 于 网 格 畸变 会 导致 的 时 间 步 缩短 ， 且 为 避免 单元 负 体积 错误 和 
网 格 异常 穿 透 需 缩小 时 间 步 长 系数 (TSSFAC ) 以 确保 稳定 性 , 实际 计算 时 长 均 高 于 SPH. 


“History Variable#6 
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8.000e-01 
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(a) Johnson-Cook 损伤 (b) Tresca 最 大 剪 应力 
图 1.6 SPH 法 模拟 990 m/s M193 冲击 6 mm Armox 600 15 us 时 刻 刘 板 破坏 情况 
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1.4 本 章 小 结 

本 章 建立 了 FEM、EFG3、EFG4、SPH 铬 芯 弹 高 速 冲击 仿真 模型 ， 分 析 四 种 单元 算 
法 的 特点 。 结 果 表 明 ，SPH 法 具有 最 佳 的 计算 效率 、 精 度 和 稳定 性 ， 能 准确 复 现 靶 板 的 
剪 切 冲 塞 破 坏 ， 因 此 以 下 仿真 均 采 用 SPH 法 进行 。 
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2 SPH 仿真 参数 优化 


2.1 SPH 方法 基本 原理 

SPH 方法 最 初 由 L. В. Глсуб ЖІК. А. Gingold65 提 出 ， 用 于 求解 天 体 物 理 问题 。 此 
方法 的 进一步 发 展 使 得 它 的 实用 性 也 扩展 到 了 固体 和 流体 力学 的 广泛 问题 上 。 为 表征 
SPH 粒子 在 时 间 步 内 的 运动 ， 需 要 构造 一 个 近似 函数 ， 其 基本 形式 为 ; 


ло) = о) Ws х) Q.D 
ИЛЛЕ WIREH ТЕУІ 
X, — X; 
W (x, R N (2.2) 


式 中 4 为 空间 维 数 ， 本 文 所 有 模型 均 为 三 维 ;，0(x) 为 辅助 了 两 数 ， 如 五 阶 样 条 函数 ， 由 参 
Ж. SPHKERN 控制 ,4 为 光滑 长 度 , 用 于 判断 附近 点 j 对 目标 点 i 的 影响 程度 。 LS-DYNA 
中 的 SPH 求解 器 使 用 可 变 的 光滑 长 度 。 通 过 各 部 件 中 节点 间距 的 最 大 、 最 小 值 和 CSLH 
常数 计算 初始 光滑 长 度 ho, HP CSLH 推荐 值 为 1.2， 以 保证 计算 域 在 各 方向 上 均 能 履 
盖 两 层 以 上 粒子 。CSLH RK, 计算 域 履 盖 面 积 越 大 。 若 值 小 于 1.05 会 导致 计算 不 稳定 ， 
大 于 1.3 会 显著 增加 计算 耗 时 。 各 粒子 具有 随时 间 变 化 的 、 独 立 的 光滑 长 度 ， 其 计算 方 
法 为 : 


а 1 
qt UOI T ROV Q.3) 


式 中 Vev 为 流动 发 散 性 ;AD 为 光滑 长 度 ， 应 介 于 最 小 值 HMINX (1-0) fl HMAX XA 
(1-0) 之 间 。HMIN 默认 值 为 0.2。HMAX 与 核 函 数 有 关 : 对 于 三 阶 样 条 函数 ， 默 认 值 
为 2.0; 对 于 五 阶 样 条 函数 , 推荐 取 3.0 或 更 大 。 以 上 参数 中 , 可 在 LS-DYNA 中 选择 的 、 
对 仿真 结果 影响 最 明显 的 是 CSLH。 此 外 还 包括 粒子 近似 理论 FORM 和 失效 粒子 控制 方 
法 了 EROD。 以 下 通过 仿真 研究 此 三 参数 对 仿真 精度 的 影响 ， 进 而 确定 最 佳 参数 组 合 。 
2.2 初始 光滑 长 度 

首先 测试 初始 光滑 长 度 参数 CSLH 的 影响 ， 网 格 /粒子 模型 及 材料 参数 与 前 文 相同 。 
分 别 对 初速 937 m/s 和 990 m/s 的 侵 彻 过 程 进行 仿真 ， 以 击 穿 / 阻 挡 情 况 、 穿 孔 和 塞 块 形 
瑶 、 全 弹 体 和 铅 蕊 所 受阻 力 评 估 其 计算 准确 性 。 参 考 推荐 值 ，CSLH 分 别 设 为 1.05、1.1、 
1.2、1.3、1.4。 两 组 CSLH=1.05 均 发 生 粒 子 超速 错误 导致 计算 异常 终止 ， 以 下 不 再 分 析 。 
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由 图 2.1 可 见 ， 弹 体 所 受阻 力 随 CSLH 增加 而 增加 ， 均 呈 增 加 一 稳定 一 衰减 趋势 ， 主 要 
差距 在 于 稳定 段 的 持续 时 间 ，CSLH 越 低 ， 稳 定 阶段 越 短 。 由 于 SPH 法 存在 拉 伸 不 稳定 
问题 , 间距 过 大 的 粒子 会 在 发 生 由 Johnson-Cook 模型 控制 的 物理 失效 前 先 断 开 连 接 ,， 即 
数值 失效 。CSLH 1.1 靶 板 背面 受 拉 伸 破 坏 较 重 ， 划 板 在 被 冲 出 塞 块 前 先 数值 失效 产生 
小 缺口 ， 而 后 剩余 粒子 受 拉 应 力 破 坏 失 效 ， 与 实验 结果 有 明显 偏差 ; 在 阻力 上 表现 为 无 
稳定 侵 彻 阶段 ， 衰 减 段 阻力 下 降 更 早 。 此 外 ，CSLH 1.1 也 唯一 在 937 m/s 时 被 击 穿 : 与 
990 m/s 类 似 , 鄞 板 背面 撕 裂 朋 落 产生 碎片 云 , 弹头 未 穿 出 , 因此 显然 是 不 准确 的 。CSLH 
1.2~1.4 数值 失效 问题 逐渐 改善 ，1.3、1.4 阻力 曲线 几乎 重合 ; 穿孔 均 为 近似 柱 形 ， 穿 出 
体 均 为 中 心 塞 块 + 环形 大 碎 块 + 弹 体 和 碎片， 与 6.5 mm Secure 500 FEM 仿真 和 实验 结果 接 
近 03]。 穿 孔 直 径 分 别 为 9.7 mm. 9.6 mm. 9.6 mm, #7 mm Armox 600 实验 结果 上 略微 偏 
高 49。。 综 上 ， 可 以 认为 ，CSLH 1.3、1.4 均 较 好 地 抑制 了 拉 伸 不 稳定 及 其 引发 的 数值 
失效 问题 ， 能 准确 地 描述 阻力 过 程 。 虽 1.4 因 计算 时 间 较 长 而 不 在 推荐 范围 内 ， 但 对 上 毁 
伤 形 貌 和 阻力 过 程 描述 更 准确 ， 以 下 仿真 均 采 用 CSLH 1.4. 


isses CSLH 1.1 Core 
сеге CSLH 1.1 Total 
CSLH 1.2 Core 
: CSLH 1.2 Total 


w 
= 
S 
1 


ға 
a 
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Resistance (kN) 


w | | 3 əн 
M 而 20 30 40 50 


Time (ps) 


(a) 937 m/s 


a 

a 

о 
1 


Resistance (kN) 
š 


Time (us) 


(b) 990 m/s 
21 不 同 初始 光滑 长 度 模拟 不 同 速度 M193 12 6 mm Armox 600 过 程 阻 力 


а 
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(a) CSLH=1.1 (b) CSLH=1.2 


23 粒子 近似 方法 
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(с) CSLH=1.3 
图 2.2 990 m/s M193 侵 彻 6 mm Armox 600 不 同 CSLH FEAR fo 


9.000e-01 
8.000e-01 
7.000e-01 
6.000e-01 
5.000e-01 
4.000e-01 
3.000e-01 


0.000e+00 | 


(d CSLH=1.4 


LS-DYNA 内 置 了 多 种 粒子 近似 方法 , 其 中 可 用 于 SMP 计算 的 方法 名 称 与 编号 如 表 
4 所 示 。Scazzosi 等 B99 的 研究 表明 FORM 4. 7. 8 不 能 求解 高 速 冲击 问题 ， 故 在 此 不 作 
讨论 。 仿 真 结 果 显 示 ， 两 种 各 向 同 性 重 整 化 近似 方法 FORM 3、6 和 两 种 各 回 异 性 近似 
方法 FORM 9. 10 JE Ж SPH-FEM 边界 处 发 生 粒 子 超速 错误 ， 导 致 计算 无 法 继续 进 


ÍT, SPH 与 FEM 无 法 通过 常规 的 固 连接 触 簿 合 ; 两 种 


兽 强 流体 方法 FORM 15. 16 均 有 


弹 体 碎片 粒子 发 生 超速 错误 ， 导致 计算 无 法 计算 进行 ,并且 相 较 于 FORM 1 SER 011 
拉 伸 不 稳定 现象 更 严重 ; 所 有 非 重 整 化 近似 方法 FORM 0、2、5、15 均 与 图 2.3 (a) 类 似 ， 
在 靶 板 SPH-FEM 边界 处 发 生 数值 失效 ，SPH 与 FEM 部 分 分 离 ， 并 有 少量 边界 粒子 和 
实体 单元 发 生 异 常 物理 失效 。 此 外 ， 从 图 2.3 (b) 终 时 刻 靶 板 毁 伤 形 貌 来 看 ， 整 体 上 ， 相 


较 于 标准 方法 ， 重 整 化 近似 方法 能 更 好 地 抑制 拉 伸 不 稳定 问题 。 


一 致 ， 近 似 方法 选择 精度 最 高 的 重 整 化 近似 FORM 1. 


4€ 2.1 LS-DYNA 4 


综 上 ， 以 下 与 前 文保 持 


可 用 于 高 速 冲击 问题 的 粒子 近似 方法 
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粒子 近似 方法 
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增强 流体 


А 
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xk 
数 


整 化 对 称 近似 
流体 近似 
重 整 化 流体 近似 


hy PRI BL 
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Ка. Ы АЕ Мыр жаа Y 
(а) FORM 0 (b FORM1 
图 2.3 重 整 化 对 弹 、 鞭 形 貌 的 影响 


2.4 失效 粒子 控制 方法 

IEROD 关键 字 用 于 控制 失效 的 SPH 粒子 。IEROD 0 表示 粒子 始终 保持 活动 状态 ， 
其 应 力 不 会 在 失效 后 自动 归 零 ， 并 且 会 保留 体积 响应 ， 被 认为 不 适用 于 大 失效 问题 ; 
IEROD 1 表示 粒子 在 失效 后 部 分 失 活 ， 应 力 自 动 归 零 并 失去 体积 响应 ， 但 仍 保留 于 计算 
域 中 以 确保 计算 稳定 性 ; IEROD 2 表示 粒子 在 失效 后 完全 失 活 ， 应 力 自动 归 零 并 失去 体 
积 响应 ; IEROD 3 保留 由 状态 方程 控制 的 体积 响应 ， 其 余 与 EROD 1 相同。 仿真 发 现 ， 
IEROD 0 与 IEROD 3 在 FEM-SPH 边界 处 计算 稳定 性 很 差 , 在 使 用 多 种 粘性 和 刚性 沙漏 
| 〈IHQ=1、2、3、4、5) 时 均 因 在 边界 处 FEM 网 格 异 常 畸 变 发 生 负 体积 错误 进而 
导致 计算 终止 。IEROD 1 与 IEROD 2 以 半 流 体 侵 御 为 主 的 阻力 上 升 和 稳定 段 几 乎 重合 ， 
仅 下 降 段 后 期 有 明显 差别 ， 最 终 塞 块 动能 分 别 为 121、211 J。 与 IEROD 1 相 比 , IEROD 
2 所 得 塞 块 主体 形状 基本 相同 , 但 以 拉 伸 破坏 为 主 的 头 部 边缘 相对 较 窄 ; 未 形成 碎片 环 ， 
而 是 一 个 较 窗 的 完整 环 与 大 量 失 效 粒 子 云 。 从 结果 来 看 ， 压 缩 破坏 过 程 对 失效 粒子 控制 
方法 不 敏感 ， 主 要 影响 在 于 拉 伸 和 剪 切 破坏 ，IEROD 2 靶 板 的 抗 拉 、 剪 韧性 相对 较 低 。 
此 外 ， 与 其 描述 相符 ，IEROD 2 计算 稳定 性 较 IEROD 1 略 差 ， 前 者 需 将 时 间 步 长 缩放 
系数 减 小 至 0.4 才能 稳定 计算 ， 否 则 会 发 生 粒 子 超速 错误 ， 而 后 者 仅 需 0.8。 综 上 ， 以 下 
仍 采 用 部 分 失 活 的 IEROD 1. 


Besse IEROD 1 


"P COEM 77 IEROD 2 
200 十 


ЕІ 
Time (ps) 


图 2.4 不 同 失效 粒子 控制 方法 模拟 990 m/s M193 {540 6 mm Armox 600 阻力 过 程 
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(а) IEROD 1 (b) IEROD 2 
图 2.5 ”失效 粒子 控制 方法 对 靶 板 毁伤 形 貌 的 影响 


2.5 Jk 


本 章 探究 了 初始 光滑 长 度 (CSLH) 、 粒 子 近 似 方法 (FORM) 和 失效 粒子 控制 方 
法 (IEROD) 对 仿真 的 影响 。 结 果 表 明 ，CSLH=1.4、FORM=1、IEROD=1 的 参数 组 合 
能 最 准确 地 模拟 铅 蕊 弹 冲 击 过 程 ， 且 具有 良好 的 计算 稳定 性 。 
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З 靶 板 硬度 对 穿 甲 过 程 的 影响 


3.1 SOARS eRe 


4340 的 性 能 与 其 热处理 方法 有 很 大 关系 , 不 同文 献 中 提出 的 数值 差别 较 大 。 由 于 文 

献 [30] 仅 给 出 了 洛 氏 硬度 (HRC)〉、 屈 服 强度 和 抗 拉 强度 ， 其 布 氏 硬度 依 GB/T 1172 换 
算 。 对 30 HRC? 302 НВУУ ЯН В5-41 侵 彻 的 研究 表明 ， 两 种 硬度 靶 板 的 抗 弹性 
能 不 存在 可 见 区 别 ,因此 本 文 所 用 4340 可 近似 为 300 HBW 的 装甲 钢 。Mars 190 属于 上 典 
型 的 大 厚度 装甲 钢 ， 便 度 随 板 厚 的 增加 而 显著 下 降 ， 本 文 及 文献 [32] 均 取 薄 板 的 硬度 上 
限 。Armox 500 和 Armox 600 应 用 较为 广泛 ， 且 均 有 明确 的 生产 标准 ， 实 测 硬 度 分 别 为 
516 НВУ 579 HBWISl。 综 上 ， 尽 管 受 限 于 生产 工艺 ， 各 靶 的 硬度 都 存在 偏差 ,但 
均 在 误差 允许 范围 内 。 
昌 于 装甲 钢 以 抗 弹 而 非 力 学 性 能 作为 指标 ， 其 强度 和 硬度 的 波动 范围 较 一 般 结构 钢 
(то, 例 。 整体 来 看 , 屈服 强度 的 数值 约 相当 于 布 氏 硬度 CBHN) 的 3 倍 。 以 下 对 300~600 
HBW 范围 内 的 5 种 钢 靶 , 分 别 在 等 靶 厚 击 穿 、 临 界 击 穿 和 半 无 限 靶 三 种 工 况 下 进行 990 
m/s M193 侵 彻 仿真 ， 以 分 析 硬 度 对 钢板 抗 铅 蕊 弹性 能 的 影响 。 表 中 列 出 了 各 材料 的 名 
称 及 标 称 性 能 。Mars 190 较为 特殊 ， 其 硬度 、 强 度 和 伸 长 率 范围 极 大 ， 并 且 可 分 为 高 低 
两 级 硬度 的 等 级 ， 分 别 记 为 Class 1 和 Class 2。 与 其 硬度 、 功 能 相近 的 Armox 370 也 分 
为 两 级 ， 均 符合 MIL-A-12560 中 对 RHA 的 规定 。 本 文 使 用 的 是 硬度 较 高 的 一 


D 


32531 ММЖФАФХЧТЕН 
材料 H/HBW omiVMPa  owMPa  ” 伸 长 率 /% 
4340 300 850 1020 14~16 
Mars 190 Class 1 240-400 600-900 800~1000 12-20 
Armox 500 480-540 1250 1450-1750 10 
Armox 600 570 - 640 1500 2000 7 


Johnson-Cook 模型 中 将 Taylor-Quinney 系数 设 为 0.9， 即 90% 的 塑性 功 转化 为 热 。 
应 变 率 在 高 速 侵 彻 中 影响 显著 ， 以 下 设 为 104/s。 图 3.1 (a) 展 示 了 等 效应 变 由 0 增长 至 2 
的 过 程 中 等 效应 变 的 变化 过 程 ， 可 见 四 种 钢 靶 的 变形 过 程 均 可 分 为 弹性 、 应 变 强化 和 热 
软化 三 个 阶段 ， 其 中 弹性 段 最 短 ， 应 变 强化 段 次 之 ， 热 软化 段 占 比 最 大 。 对 最 大 等 效应 
力 与 硬度 的 关系 进行 线性 拟 合 ， 结 果 如 式 (3.1) 所 示 ， 误 差 如 图 3.1 (Ы). 
=3.528H +80.54 (3.1) 


бос тах 
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= m N 
CD оо N 
сл о л 
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Stress (МРа) 
o 
S 


(a) 应力 一 应 变 关系 


A 
е 
о 
б 


m Simulation 


Linear 


Maximum stress (MPa) 


Hardness (HBW) 
(b) ”最 大 应 力 一 硬度 
图 3.1 考虑 应 变 率 强化 和 热 软 化 效应 的 应 力 一 应 变 关系 


钢 的 硬度 与 强度 的 关系 有 多 种 经 验算 法 。Tate 等 多 认为 在 长 杆 高 速 侵 彻 中 靶 板 材料 
可 近似 为 理想 弹 塑性 ， 其 屈服 强度 oy.Tate 相当 于 4.2 倍 布 氏 硬 度 。JC 模型 线性 拟 合 的 截 
距 仅 相当 于 oycmax 的 不 足 10%, 在 粗略 计算 时 可 近似 为 正比 ,但 Tate 式 的 斜率 明显 偏 高 ， 
以 致 过 度 高 估 了 装甲 钢 的 强度 。Rosenberg 等 以 准 静 态 一 常温 条 件 下 e=0.27 时 的 应 力作 
为 等 效 强 度 沪 ]， 由 于 同样 久 视 了 热 软 化 ，oRes 在 数值 上 与 оутае 接近 。 偏 差 可 能 来 自 于 
两 者 研究 对 象 的 不 同 : Tate 等 和 Rosenberg 等 所 用 邯 涵 盖 铝 、 铀 、 钢 等 各 种 金属 ， 硬 度 
普遍 在 300 HBW 以 下 ; ПЖ 262488, Н.Н, > 300HBW。 此 外 ， 各 金属 强度 
标准 均 存 在 硬度 一 抗 拉 强 度 插 值 换算 方法 ， 以 GB/T 1172 为 例 C9， 此 换算 方法 除 基 础 的 
硬度 外 还 考虑 了 钢材 化 学 成 分 的 影响 。 以 铬 镍 钢 估算 300 HBW 的 4340， 以 超 高 强度 钢 
估算 400-600 HBW 1, SERU om 所 示 。 尽 管 抗 拉 强 度 测 试 中 应 力 随 应 变 增 加 而 下 
降 的 原因 是 颈 缩 ， 而 与 热 软 化 无 关 ， 并 且 忽 略 了 应 变 率 强 化 效应 ， 但 两 者 呈现 出 良好 的 
相关 性 ， 各 钢 训 的 偏差 均 在 10% 以 内 。 
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3532 等 效 强度 估算 方法 比较 
等 效 强度 /MPa RAIN 
OjCmax OyTate ORos Om Oft Om 
300 1152 1260 1198 1046 1.1 92 
400 1450 1680 1448 1370 2.9 5.5 
500 1888 2100 2101 1716 23 9.1 
600 2182 2520 2512 2180 07 041 


分 别 以 400 和 500 HBW 级 装甲 钢 作为 验证 :Hardox 400 的 orcmax 为 1484 MPa"; 
另 一 Armox 500 模型 为 1929 МРа/44, 30 PM 为 1993 MPal!, Protac 500 为 1991 MPa", 
考虑 到 装甲 钢 的 硬度 和 强度 普遍 存在 一 定 生产 误差， 如 Armox 560Т 标 称 最 小 屈服 强度 
1300 MPa®!!, 51] 1711 МРа2!; Demox Protect 500 标 称 屈服 强度 1500 MPaP?l, 51] 
1592 MPa'?l; Hardox 400 标 称 屈服 强度 1100 МРа/54, “2 | 1148 MPal?!, 经验 式 的 偏差 
在 合理 范围 内 。 综 上 ， 对 于 硬度 在 300~600 HBW 范围 内 的 装甲 钢 ， 可 由 Johnson-Cook 
模型 或 式 (3.0D) 估 算 其 在 侵 彻 过 程 中 的 最 大 应 力 ， 也 可 由 式 G3.D) 反 向 初步 评价 
Johnson-Cook 模型 的 可 靠 性 。 

为 研究 硬度 的 具体 影响 , 以 下 对 等 厚度 击 穿 、 弹 道 极 限 、 半 无 限 靶 三 种 工 况 下 M193 
的 侵 彻 进行 测试 ， 初 速 均 为 990 m/s. 


3.2 ŞEREF 


为 确保 各 组 均 能 被 击 穿 ， 取 和 靶 板 厚度 为 600 HBW 的 临界 击 穿 厚度 6 mm. Hë 3.3 
可 见 ， 随 着 硬度 的 增加 ， 靶 板 的 破坏 形式 由 拉 伸 逐渐 过 渡 至 剪 切 ， 穿 透 后 弹 体 前 的 靶 板 
剩余 厚度 逐渐 增加 。300 和 400HBW 表现 为 拉 伸 破坏 ， 靶 后 雁 片 各 方向 分 布 较为 均匀 ， 
整体 上 附着 于 残余 弹 体 上 ， 旦 近似 半 椭 球状 ，500 和 600HBW 表现 为 剪 切 破坏 ， 其 中 
500 HBW 由 于 韧性 较 好 ， 冲 出 的 塞 块 较 600 HBW 时 完整 。 

依照 发 生 部 位 ， 靶 板 的 破坏 可 分 为 正面 压缩 、 环 形 前 切 和 背面 拉 伸 三 部 分 : 

正面 破坏 来 源 于 弹 对 部 的 直接 挤 压 ， 四 强度 靶 板 正面 破坏 均 仅 限 于 很 薄 的 一 层 ， 表 
明 在 6 mm 时 抗 压 不 是 抗 弹性 能 的 主要 因素 。 正 面 破坏 可 分 为 中 心 和 边缘 两 部 分 ;中心 
损伤 区 随便 度 的 提高 而 厚度 减 小 ; 边缘 损伤 区 随便 度 的 提高 而 增 大 。 尽 管 发 生 完全 破坏 
的 范围 很 小 ， 如 靶 体 损伤 一 栏 所 示 ， 但 类 似 于 长 杆 高 速 侵 彻 ， 塑 性 流动 也 会 造成 抗 阐 能 
力 的 损失 。 由 于 半 无 限 训 高 度 承 压 ， 其 应 力 三 轴 度 很 低 ， 靶 板 的 失效 仪 局 限于 轴线 处 ， 
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但 因 纯 粹 塑性 而 开 孔 。 因 此 , 以 侵 彻 结束 时 的 残余 靶 厚 评估 实际 正面 破坏 程度 , 其 中 300 
和 400 HBW 由 于 背面 破坏 较 强 

环形 破坏 类 似 于 典型 的 剪 切 种 塞 , 但 因 误 板 存 在 变形 而 附加 了 拉 伸 的 影响 。 由 表 3.3 
中 部 板 损伤 和 失效 两 行 可 见 ， 典 型 的 环形 破坏 面 仅 存 在 于 400-600 НВУ. В тя 
度 的 减 小 ， 剪 切 破坏 由 600 HBW 时 的 近似 柱 面 逐渐 扩张 为 400 HBW 的 约 45” 锥 面 ， 
同时 失效 面 的 厚度 逐渐 减 小 。300 HBW 由 于 受 正 面 破坏 过 多 ， 误 板 向 背面 显著 弯曲 ， 
因此 受 虽 没有 剪 切面 ， 但 由 临 穿 透 损伤 可 观察 到 拉 伸 破坏 。400 HBW 也 存在 相似 现象 。 

背面 破坏 来 源 于 凸 出 、 弯 曲 导 致 的 拉 伸 ， 分 为 表面 和 轴线 两 部 分 。 轴 线 破 坏 来 源 于 
弹 尖 的 集中 拉 伸 ， 表 面 破坏 来 源 于 整体 弯曲 。 刘 板 硬度 和 强度 越 高 ， 抗 弯 性 能 越 好 ， 同 
时 残余 靶 厚 越 大 ， 整 体 刚 度 越 高 背面 变形 越 小 ， 也 就 越 不 容易 发 生 拉 伸 破坏 。 也 与 高 应 
力 三 轴 度 下 的 塑性 有 关 。 例 如 对 于 400 HBW， 其 本 身 抗 弯 性 能 不 佳 ， 且 高 应 力 三 轴 度 
下 塑性 很 低 ， 因 此 背面 破坏 最 为 严重 ，300 HBW 虽然 抗 弯 强度 不 及 400 HBW， 但 塑性 
更 好 ， 因 此 表面 破坏 较 轻 。 


X33 Яш 
硬度 /HBW 300 400 500 600 


穿孔 形 貌 


临 穿 透 损伤 
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B FUT ° 
Е : : 
由 表 3.4 可 见 , 穿孔 直径 整体 上 与 靶 板 硬度 负 相 关 ， 主 要 原因 是 随 着 靶 板 硬度 提高 ， 
铅 亿 所 受 轴 疝 阻力 也 相应 增强 ， 进 而 铅 、 铜 径 向 流动 导致 的 扩张 越 强 。 


600 HBW 与 此 


趋势 异常 是 因为 靶 板 硬度 足够 高 ， 以 至 于 在 发 生 剪 切 冲 塞 前 几乎 不 受 侵蚀 ， 并 仅 有 轻微 


变形 。 


因此 最 终 的 穿孔 直径 更 接近 于 冲 塞 直径 ， 而 不 是 其 它 靶 的 侵蚀 直径 。 


残余 靶 厚 与 


辑 板 硬度 正 相 关 ， 即 靶 受 侵 s MON 与 硬度 负 相 关 。 其 中 300 HBW 虽然 弹 体 前 


方 没有 剩余 肢体 ， 但 部 并 非 被 完全 侵 


由 表 3.3 中 穿孔 形 貌 一 栏 可 见 ， 因 为 残余 弹 体 
的 动能 较 高 且 头 部 较为 尖锐 ， — E 拉 断 。 两 侧 的 残余 靶 体 虽然 


ENS 


厚度 不 均匀 ， 但 从 图 像 可 定性 判断 薄 于 400 HBW， 即 与 前 文 所 述 趋势 相同 。 


253.44 990 m/s M193 临界 击 穿 数 值 
人 硬度 /HBW 穿孔 直径 /mm ”残余 靶 厚 /mm 动 


300 8.88 
400 9.03 
500 9.53 
600 8.9 


3.3 弹道 极限 


由 表 3.4 AL, Bae 
的 环形 破坏 面 。 破 坏 面 由 


300 HBW 的 约 45° 
表明 驾 板 在 冲 塞 前 发 生 的 弯 拉 变形 逐 


0 
1.03 
4.00 
4.60 


结果 

能 消耗 订 
1353 
1594 
1707 
1662 
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板 在 临界 击 穿 时 均 发 生 剪 切 冲 塞 ， 击 穿 损 伤 一 栏 中 均 可 见 明 下 
截 锥 形 逐 渐 收 缩 为 600 HBW 的 近似 柱 形 ， 
簿 渐 减 弱 。300 HBW 在 冲 塞 的 同时 发 生 与 6 mm 相似 


== 
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的 拉 伸 破坏 ， 塞 块 与 靶 板 出 口 处 的 翻 展 被 拉 断 形成 大 碎 块 ; 400 HBW SEXE RINL7J —5 
度 下 塑性 较 低 ， 发 生 背 面 朋 落 并 产生 完全 失效 的 碎 导 ， 以 致 塞 块 较 薄 ; 500 和 600 HBW 
为 典型 冲 塞 ， 其 中 500 HBW 由 于 塑性 较 高 ， 塞 块 更 为 完整 。 各 计 临 界 击 穿 厚度 如 表 3.6 
中 Ti penmax 所 示 。 


323.55 M Fed л РЕЗ 
300 HBW 400 HBW 500 HBW 600 HBW 


Tn. 


SERE AN 


PEIR Ø 


击 穿 破 坏 


% 


ж 
П 


м 
П 


A 
Й 


т 


Maximum penetration thickness (mm) 


300 400 500 600 
Hardness (HBW) 


32 RATATE 


Ds 
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由 表 3.6 可 见 ， 背 面 变 形 量 与 残余 邯 厚 与 台 板 硬度 分 别 呈 负 相关 和 正 相 关 ， 表 明 随 
着 硬度 的 提高 ， 主 要 破坏 形式 由 弯曲 一 拉 伸 转变 为 剪 切 。400 HBW 由 于 发 生 骨 落 ， 残 
余 厚 度 较 低 。 各 误 板 破坏 形式 均 与 临界 击 穿 相似 ， 均 存在 环形 破坏 面 : 300 HBW 面 厚 
度 最 大 ， 背 表面 存在 损伤 但 未 完全 失效 ，400 HBW REMI, DUE TUAE 
常 偏 低 ;500 和 600 HBW 靶 损 伤 形式 相同 ， 均 包括 环 状 剪 切 和 背面 弯曲 两 部 分 。 


表 3.5 ЯН Сл PSERUEA 


300 HBW 400 HBW | 500 HBW | 600 HBW 
| | š ІШІ TA L : ІШІШІШ/ 
损伤 
EE 
损伤 (复原 ) | E 
| i 
#36 背面 最 大 变形 量 
 HJWBW Гош BFD/mm Гиш. 
30 | 9 57 471 
400 125 5.55 3. 86 
500 6.5 5.17 5,13 
600 6 4, 94 5.10 
3.4 ІП 


为 明确 半 流 体 侵 彻 性 能 ， 进 行 半 无 限 鞭 仿真 。 取 靶 厚 25 mm 以 完全 排除 背面 变形 、 
应 力 波 反 射 等 背面 效应 的 影响 。 由 图 3.3 可 见 ，300 和 400 HBW 靶 均 受到 显著 侵蚀 ， 弹 
体 雄 片 形成 明显 管状 残余 体 ， 靶 板 形成 翻 唇 ，500 靶 翻 唇 在 终 时 刻 已 断裂 ， 向 后 随 弹 体 
碎片 一 同 飞 散 ， 形 成 靶 前 碎片 云 ，600 HBW 靶 仅 受 轻 微 侵蚀 。 如 图 3.5 Bray, Ara SOS 
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发 生 球面 和 轴线 两 部 分 失效 : 300. 400 和 500 HBW 车 的 球面 失效 区 域 均 深 入 靶 体 内 部 ， 
在 其 外 部 存在 无 失效 区 ， 即 图 3.3 与 3.4 中 的 翻 唇 /大 碎片 ，600 HBW 由 于 硬度 高 且 塑 性 
低 ， 球 面 失 效 厚 而 浅 ， 不 存在 外 部 无 失效 区 ， 因 此 未 能 形成 翻 唇 。 


(a) 300 HBW (b) 400 HB (с) 500 НВУ (а) 600 HBW 
图 3.3 终止 时 刻 弹 靶 形 貌 


| | Тт; | 六 Ш.Д [M 
(a) 1.75 us (b) 2.24 us (c) 2.72 us (d) 3.12 us 
图 3.4 500 HBW Н/Е. Wr ES SE BUE FT KAOLE 
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(a) 300 HBW (b) 400 HBW (с) 500 HBW (d) 600 HBW 


图 3.5 终止 时 刻 靶 体 损 伤 


由 头 尾 速度 可 分 为 三 个 阶段 ， 以 下 以 侵 彻 深度 最 大 的 300 HBW 为 例 : 

1) 初始 瞬 态 阶段 : 约 0~0.5 mm。 此 阶段 弹 鞭 开始 接触 。 由 于 弹 体 强度 远 低 于 靶 板 ， 
接触 面 上 的 弹 体 迅 速 变形 、 扩 张 ， 头 部 速度 迅速 下 降 

2) 准 定常 阶段 ， 对 应 侵 彻 深度 约 0.5-4.2 mm。 此 阶段 类 似 于 半 流 体 长 杆 侵 彻 的 主 
112 2. 尾 速度 及 两 者 差 值 近似 保持 恒定 ， 其 中 头 部 速度 约 在 100 
m/s 左右 ， 尾 部 速度 接近 于 其 初始 速度 990 m/s。 根 据 实 际 发 生 侵 彻 的 部 件 ， 可 分 为 两 个 
小 阶段 : 


2.1 ”被 甲 侵蚀 阶段 ， 0.5~1.3 mm。 此 阶段 弹 体 强 度 较 高 的 被 甲 部 分 与 靶 体 接 
fih 1.3 mm 时 尖端 被 甲 完 全 破坏 ， 头 部 速度 出 现 谷 值 ， 头 尾 速度 差 出 现 峰 值 ， 

22 ” 铝 蕊 侵蚀 阶段 ，1.3~4.2 mm。 此 阶段 为 铅 芯 弹 侵 彻 半 无 限 靶 的 最 主要 阶 
段 ， 铅 芯 在 此 阶段 与 间 板 在 轴 向 直接 接触 ， 被 甲 向 外 翻 折 ， 仅 对 铅 芯 有 微弱 的 径 向 
约束 作用 。 

3) 停止 阶段 : 4.2~5.5 mm， 相 当 于 半 流 体 长 杆 侵 彻 的 侵 彻 末期 。 此 阶段 弹 体 不 再 满 
足 流体 头 部 一 刚体 主体 假设 。。 可 分 为 三 个 小 阶段 。 

3.1 42-5.0 mm: 头 部 速度 仍 保持 在 约 100 m/s 不 变 ， 尾 部 速度 开始 下 降 。 
同时 靶 体 扩张 减缓 ， 至 5.0 mm 时 弹 体 直径 达到 最 大 值 ， 尾 部 减速 度 逐 渐 升 高 

3.2 50-5. mm: 尾部 快 减速 ， 阻 力 下 降 。 至 5.5 mm Sk ABR REVS, 30 
体 不 再 受 破坏 。 

33 5.5~ 终 止 ; 此 阶段 为 靶 体 的 弹性 回 弹 ， 弹 坑 深 度 略 微 减 小 ， 侵 彻 深度 随 
之 振荡 。 此 阶段 相对 于 整个 侵 彻 过 程 而 言 占 比 极 小 。 
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一 一 300HBW Head 
一 一 300HBW Тай 
— — A00HBW Head 
一 一 400HBW Тай 
一 一 500HBW Head 
MEMO тама M 一 一 500HBW Tail 
一 一 600HBW Head 
一 一 600HBW Tail 


1000 -~ 


Velocity Difference (m/s) 
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(a) 头 尾 速度 


) 
N 
3 
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N 
s 
1 


Depth of Penetration (mm) 


(b) 头 尾 速度 差 


Depth of Penetration (mm) 


(c) 阻力 
图 3.4. M193 侵 彻 各 硬度 钢 靶 过程 
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(a) DoP=0.5 mm (b) DoP=1.3 mm (с) DoP=4.2 mm(d) БоР-5.0 mm 
图 3.5 300 HBW 弹 靶 变形 与 破坏 


尽管 各 钢 靶 均 可 观察 到 半 流 体 侵 彻 行为 ， 但 由 于 铬 必 弹 长 度 较 短 ， 且 不 由 单一 材料 
组 成 ， 其 侵 彻 过 程 中 速度 和 阻力 的 变化 与 典型 长 杆 弹 有 较 大 差别 。 由 图 3.4 可 见 ， 头 部 
速度 存在 较 大 幅度 的 波动 ， 整 体 上 与 靶 体 便 度 呈 负 相关 ， 同 时 头 尾 速度 差 与 间 体 硬度 正 
相关 ， 表 明 提 高 轩 板 硬度 能 更 有 效 地 损耗 弹 体 ; 阻力 在 整个 过 程 中 不 存在 稳定 阶段 ， 而 
是 始终 随 弹 体 的 钝 化 和 扩张 而 增 大 。 

由 图 3.5 可 见 ， 侵 彻 深度 与 最 大 阻力 均 与 靶 体 便 度 呈 近 似 线性 关系 。 与 3.1 节 对 最 
大 应 力 的 分 析 相 同 ，400 和 500 HBW 撮 的 抗 弹 能 力 分 别 较 拟 合 线 俩 低 和 侦 高 ， 表 明寺 
体 的 材料 模型 与 硬度 存在 误差 。 
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240 ! 1 Т 1 
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=, 
(b) 最 大 阻力 


图 3.6 990 m/s M193 侵 彻 300~600 HBW ЖИН £i HL 


3.5 等 效 系 数 

分 别 以 对 半 无 限 靶 侵 彻 深度 (Kr) 、 对 有 限 厚 靶 的 临界 穿 透 厚度 (Kp) 衡量 相对 于 
300 HBW 车 的 抗 弹性 能 ， 并 最 大 应 力 (Ko) 辅助 比较 。Borvik 等 对 尖 头 刚体 弹 的 实验 
结果 表明 ， 苯 板 的 抗 弹性 能 仅 与 有 关 ， 并 提出 装甲 钢 的 强度 远 比 塑性 重要 01。 然 而 ， 由 
K 3.7 和 图 3.7 可 见 , 三 种 等 效 方 法 所 得 结果 有 很 大 偏差 : Kr 与 Ke 均 呈 现 与 硬度 的 良好 
相关 性 ， 表 明 对 于 仅 少 量 失 效 的 半 无 限 襄 抗 弹性 能 与 训 板 硬度 /强度 高 度 相 关 ; Kp 尽管 
同样 与 硬度 正 相 关 ， 但 400 HBW 明显 较 其 余 偏 低 ， 而 300. 500 和 600 HBW 5 5; 55 9 
等 效 值 相似 的 线性 关系 。Kr 与 Ke 的 偏差 主要 来 源 于 鞭 板 破坏 形式 的 差异 ， 对 于 临界 击 
穿 的 有 限 厚 靶 ， 其 塑性 决定 了 抗 拉 伸 和 剪 切 的 能 力 ， 其 中 对 于 强度 较 低 、 弯 曲 幅 度 较 大 
的 300 和 400 HBW 以 前 者 为 主 ， 对 于 强度 较 高 、 弯 曲 幅 度 较 小 的 500 和 600 HBW 以 后 
4A . 400 HBW 单 轴 拉 伸 失 效应 变 1.38, 双 轴 拉 伸 失效 应 变 0.40; 作为 对 比 , 300 HBW 
单 轴 拉 伸 失 效应 变 3.35， 双 轴 拉 伸 失 效应 变 2.93. 400 HBW 极 低 的 双 轴 拉 伸 失效 应 变 
导致 其 更 易 破坏 ， 因 此 需要 较 预 期 更 大 的 厚度 以 抵消 此 缺陷 的 负面 影响 。 


表 3.7 三 种 等 效 方法 比较 
H/HBW Kr Ke K, 
300 1 c3 l 
40 12 128 126 
500 1.38 2.07 1.64 
600 1.5 2.60 1.89 
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—m— strength 
—e— thickness 


Effectiveness 
A 
N 
o 


300 400 500 60 
硬度 (HBW) 


图 3.7 三 种 等 效 方法 比较 


3. 6 本 章 小 结 

本 章 基于 2 节 优化 的 SPH-FEM 耦合 模型 ， 比 较 了 等 厚度 击 穿 、 弹 道 极限 和 半 无 限 
靶 三 种 训 体 厚度 情况 下 铅 芯 弹 的 侵蚀 行为 。 结 果 表 明 : 铝 芯 弹 对 衣 板 的 破坏 包含 半 流 体 
侵 彻 和 前 切 冲 塞 两 部 分 ， 其 中 中 低 硬度 肢 以 前 者 为 主 ， 高 硬度 靶 以 后 者 为 主 。 但 相 较 于 
典型 半 流 体 侵 彻 ， 铅 芯 弹 冲 击 存在 非 均 质 弹 体 、 背 面 效应 等 非常 规 因 素 ， 因 此 侵 彻 过 程 
中 的 速度 、 阻 力 等 量 的 变化 存在 差别 ， 而 相 较 于 典型 冲 塞 ， 铅 芯 弹 冲击 存在 弯曲 变形 、 
拉 伸 断 裂 等 异常 现象 ， 以 致 丢 板 的 毁伤 形 貌 与 理想 情况 不 同 。 


-28- 


EE ТА A ҒА, 


4 理论 建 模 


综合 考虑 低 强 度 革 的 半 流 体 侵 彻 和 高 强度 对 的 冲 塞 两 种 破坏 形式 ， 本 文采 用 半 流 体 
侵 和 初 体现 弹 靶 相互 侵蚀 的 过 程 ， 以 剪 切 冲 寒 判 定 击 穿 。 由 于 低 硬度 肢 的 破坏 除 剪 切 冲 蹇 
外 也 包括 拉 伸 ， 而 根据 Johnson-Cook 模型 ， 靶 板 材料 的 抗 拉 卦 性 显著 弱 于 抗 剪 ， 故 此 模 
型 求 得 的 对 低 硬度 靶 的 极限 穿 透 厚度 会 较 实际 值 偏 低 。 此 外 ， 馈 蕊 步枪 绊 几 何 外 形 与 理 
想 钝 头 长 杆 有 较 大 区 别 ， 且 密度 在 轴 癌 和 径 向 上 分 布 缘 不 均匀 ， 不 能 使 用 简单 的 表达 式 
精确 描述 ， 因 此 本 模型 仅 用 于 定性 解释 和 分 析 侵 彻 过 程 。 
4.1 ЖИ 

Alekseevskii-Tate 半 流 体 侵 彻 模 型 提出 于 1967 5541, 首次 将 固体 的 高 速 侵 彻 与 流 
体 动力 学 相 结合 ， 是 长 杆 高 速 侵 彻 最 为 经 典 的 理论 模型 ， 也 是 后 世 绝 大 多 数理 论 模 型 的 
基础 。 由 于 A-T 模型 仅 体现 一 维 效应 ， 且 对 侵 彻 过 程 的 描述 过 于 平均 化 ， 其 计算 精度 较 
低 。 国 内 外 学 者 对 A-T 模型 进行 了 了 诸多 改进 ， 其 中 最 具有 代表 性 的 是 Rosenberg 和 
Walker-Anderson 模型 。Rosenberg 等 注意 到 侵 彻 过 程 中 弹 靶 界面 压力 从 弹 体 头 部 中 心 辐 
边缘 的 衰减 S11。 为 体现 此 二 维 效 应 ， 在 А-Т 模型 的 基础 上 引入 刚性 杆 体 和 和 钝 化 头 部 / 诅 
体 等 效 截面 积 的 概念 。Rosenberg 模型 的 基本 形式 为 : 


1 1 
Л -v,y +Y, -S PMs +R) 


dv, 1, 

» pl (4.1) 
"um T WE M; 

dt 

en 

dt 


式 中 pp I pi 2) 9] 73 CRUEL BE уы vi 分 别 为 弹 体 头 部 和 尾部 速度 ; Y, ОЙЫ Kun S, 
Ri ARRIETA; S 为 弹 、 误 等 效 截面 积 之 比 ; 1 为 弹 体 剩余 长 度 ，p 为 侵 彻 深度 。 式 中 
Y». R. Si 需 由 经 验 估算 或 从 仿真 结果 拟 合 。Si 参 考 Rosenberg 等 的 理论 分 析 中 1 和 王猛 等 
B59 的 实验 结果 ， 对 于 一 般 长 杆 弹 约 为 2; TE A-T 及 其 它 不 考虑 二 维 效应 的 模型 中 为 1。 
四 个 方程 中 最 主要 的 为 首 个 方程 ， 通 常 也 被 称 为 平衡 方程 ， 用 于 计算 va: 后 三 个 方程 分 
别 计算 弹 体 减 速 、 弹 体 侵蚀 各 体 侵蚀 。 由 于 铅 蕊 和 被 甲 的 强度 均 远 低 于 四 种 装甲 钢 ， 
对 的 取 值 不 会 对 模型 求解 造成 显赫 影响 ， 初 设 为 50 MPa。 因 此 ， 模 型 中 的 未 知 量 为 等 
效 截 面积 比 % AREA EH 7] Ri。 
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及 可 通过 界面 击溃 现象 求解 ， 即 由 w=0 时 的 平衡 方程 一 一 式 (4.2) 一 一 求解 S 和 Ri。 
3.3 和 3.4 177507389, 300-500 HBW 和 600 HBW 部 的 损伤 形式 存在 明显 差异 ， 其 中 后 
者 在 临界 击 穿 时 塞 块 几乎 无 弯曲 ， 在 半 无 限 靶 时 不 产生 翻 恨 。 因 此 ， 当 靶 板 硬度 在 
500-600 HBW 之 间 的 某 一 数值 时 发 生 界 面 击溃 ， 即 侵蚀 速度 的 下 限 高 于 弹 体 初 速 990 
m/s. 


l 
22-58, (4.2) 


Rosenberg 模型 中 Ri 由 式 (4.3) 计 算 。 式 中 01 为 靶 体 等 效 强度 ， 对 于 500 HBW 417) 
1637 MPa， 对 于 600 HBW 靶 为 1799 MPa; Ei 为 对 体 弹性 模 量 ， 对 于 钢 为 206 GPa; v 
为 泊 松 比 , 对 于 钢 为 0.3。 由 以 上 参数 计算 , WJ 500 和 600 HBW 的 Re 3073 4772 ЖІ 5146 
MPa。 考 虑 到 500 HBW 15006091, Н 600 HBW 靶 表 面 仍 有 少量 失效 ， 初 取 临 界 界 
面 击溃 时 R=5000 MPa, Jl] S=1.11， 显 著 小 于 钨 合金 长 杆 弹 侵 彻 钢 靶 的 等 效 截 面积 ， 表 
明 由 于 铬 芯 弹 速度 慢 、 强 度 极 低 ， 其 径 向 流动 、 侵 蚀 能 力 很 弱 。 

с, VBE, 

Bo Ane 
ЖЕЗ) PSY ЛИЕ УРЕ Е ДАН НЕН BKK, CATT, RE ER ITE 

圆柱 的 铜 一 铅 结构 近似 为 全 铅 圆锥 一 圆柱 组 合体 。 对 于 M193 8: 

076 0<1<7.5 

| 57 75<1<193 

以 上 式 计算 990 m/s M193 侵 彻 半 无 限 靶 过 程 ，300~500 HBW 驾 侵 彻 深度 分 别 为 
3.79. 2.84. 0.86 mm， 准 定常 段 平 均 头 部 速度 vn 4) 807J 162.2. 127.1. 42.6 m/s. BH 
际 侵 彻 半 无 限 靶 过 程 相 比 ， 理 论 计 算 所 得 侵 彻 深度 明显 偏 低 ; 平均 vs 由 于 仿真 结果 存在 
大 幅 波 动 ， 难 以 定量 判断 ， 但 其 范围 基本 准确 。 偏 差 主要 来 源 于 弹 体 结构 的 简化 和 
Rosenberg 模型 自身 缺陷 。 受 限于 长 杆 半 流 体 侵 彻 研究 的 局 限 性 ， 目 前 尚 无 考虑 杆 体 非 
理想 结构 的 理论 计算 模型 。 尽 管 本 文中 铅 蕊 和 被 甲 的 强度 尽管 均 远 低 于 靶 钢 ,但 两 者 间 
仍 存 在 很 大 差距 ， 因 此 平均 杆 体 强度 的 选取 依赖 于 经 验 ， 不 能 通过 理论 计算 直接 获得 。 
张 清 源 的 对 照 性 研究 表明 ，Rosenberg 模型 本 身 预 测 侵 初 深度 偏 低 S33。 两 者 的 共同 偏差 
导致 半 流体 侵 彻 过 程 无 法 准确 计算 ， 但 侵 彻 深度 来 看 能 定性 反映 靶 板 的 抗 弹 能 


4.2 Зз 


对 于 驾 体 几乎 不 受 侵蚀 ， 以 及 穿 甲 末期 靶 板 背面 显著 凸 出 ， 不 再 可 被 视 为 半 无 限 刘 
的 情况 ， 采 用 Recht 动量 一 能 量 平 衡 理论 判定 是 否 发 生 冲 赛 穿 孔 。Recht 理论 的 基本 假 


R = ] (4.3) 


(4.4) 


-20- 


EE ТА A ҒА, 


设 为 : 1) 弹 体 和 靶 板 材料 不 同 ， 因 而 声波 波 速 不 同 ，2) 弹 体 不 是 刚性 的 ; 3) 弹 体 和 
塞 块 的 径 向 尺寸 可 能 不 同 。 高 速 铅 蕊 弹 冲 击 钢 贰 与 以 上 假设 均 相 符 。Recht WA, 136 
穿孔 过 程 实际 上 是 一 个 动能 转换 的 过 程 ， 弹 体 原 有 动能 Exo 转化 为 : 

1) 弹 体 剩余 动能 和 塞 块 动能 ， 即 前 文 的 穿 出 体 动能 Ex: 

2) 弹 孔 周围 靶 体 材料 剪 切 破坏 吸收 的 应 变 能 Ws 

3) 弹 驾 伴 撞 过 程 中 形成 共同 速度 消耗 的 能 量 Wye 
由 此 得 到 能 量 守恒 方程 : 


E, = E, +W,+W, (4.5) 
由 于 弹 体质 量 损失 和 新 体 的 压缩 变形 和 侵蚀 的 主要 过 程 已 在 半 流 体 侵 彻 部 分 计算 ， 
且 剪 切 破 坏 耗 时 相对 较 短 ， 故 假设 冲 塞 前 后 弹 体 质量 不 变 ， 塞 块 质量 由 剩余 靶 板 厚度 和 
弹 体 直径 确定 ; 球 仅 包括 剪 切 失效 应 变 能 。 因 剖 塞 发 生 时 弹 体 尖端 已 被 完全 匆 粗 ， 故 近 
似 弹 体 为 均匀 圆柱 体 。 由 此 得 到 实际 冲 塞 过 程 的 能 量 守恒 方程 : 

1 

D 
AP: mp、mi 分 别 为 弹 体 、 塞 块 质 量 ， 可 由 半 流 体 侵 彻 模型 粗略 确定 ，Di 为 蹇 块 直径 ， 
参考 仿真 与 半 流 体 理论 分 析 结 果 粗 略 计 算 时 取 弹 体 扩张 后 直径 , BY 1.5 倍 弹 体 初始 直径 ; 
vos vr АНАЈА. BIRERE, HE vo 由 半 流 体 侵 彻 模型 确定 ， 对 于 杆 体 〔( 铅 蕊 〉 强度 
xf T REDI M193 可 取 初 速 作为 近似 ; ti 为 肢体 残余 厚度 ， 由 半 流 体 侵 彻 模型 确定 ; 
六 为 剪 切 应 变 区 径 向 等 效 厚度 ， 与 靶 体 塑 忻 有 关 ， 需 由 仿真 结果 反 辐 推算 ;ws 为 剪 切 失 
效应 变 能 密度 ， 由 Johnson-Cook 材料 模型 确定 。 辱 以 上 方程 有 对 六 的 正 实数 解 ， 则 判 
XE REA AM H E o 
以 990 m/s M193 临界 击 穿 6 mm Armox 600 靶 工 况 验 证 以 上 剪 切 冲 塞 模型 。 为 确保 
准确 性 ， 冲 塞 破 坏 初 步 以 更 精确 的 仿真 结果 计算 。 以 弹头 双 圆 驳 部 被 完全 侵 彻 作为 冲 塞 
判定 起 始 。 由 于 弹 体 尾 部 速度 在 半 流 体 侵 彻 过 程 几 乎 无 下 降 , 假设 弹 体 速度 wo 在 冲 塞 开 
台 时 仍 为 990 m/s. 其 余 参 数 分 别 为 : 塞 块 与 弹 体 的 整体 速度 为 186 m/s, mp =2.4 g, m=2.2 
о, D=8.6mm, t=4.8 mm， 解 得 r=1.0 mm。 尺 管 Recht 模型 理论 上 能 较为 准确 地 判定 
剪 切 冲 塞 ， 但 弹 体质 量 mp 的 计算 依赖 于 较 精 确 的 几何 模型 。 知 选用 半 流 体 侵 彻 模 型 中 
简化 的 全 铅 柱 体 ， 则 上 述 工 况 中 的 六 将 增 大 至 2.6 mm， 较 实际 情况 显著 偏 高 。 


ul j mm , 
v =—(m, +m,)v. - z Dro, +—— Vo (4.6) 
DER 4 ` ` 2m +m, 
p 
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43 ЖЕ МА 

本 章 基 于 第 3 节 仿真 结果 ， 建 立 了 简化 的 Rosenberg 半 流 体 侵 彻 一 Recht 剪 切 冲 塞 
模型 。 受 限于 铅 芯 一 铜 被 弹 复 杂 的 结构 和 理论 模型 自身 的 简化 与 不 准确 性 ， 尽 管 近似 理 
论 模型 可 定性 解释 其 穿 甲 原理 ， 但 精确 判断 击 穿 仍 依 赖 于 仿真 。 


—= 
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本 文通 过 FEM、EFG、SPH 三 种 方法 研究 了 高 速 铅 忌 步枪 弹 侵 彻 不 同 强度 钢 刘 的 行 
为 ， 探 索 了 初始 光滑 长 度 参 数 、 粒 子 近似 方法 和 失效 粒子 控制 方法 对 SPH-FEM 耦合 方 
法 仿真 精度 的 影响 ,基于 半 流 体 侵 彻 与 剪 切 冲 塞 模型 总 结 了 高 速 铅 芯 步枪 弹 侵 彻 钢 靶 的 
规律 ， 主 要 结论 如 下 : 

(1) 以 SPH 模拟 发 生 大 变形 和 失效 的 靶 板 核心 区 与 弹头 , 以 FEM 法 模拟 仅 小 变形 、 
受 固 定 约束 的 靶 板 外 周 区 的 SPH-FEM 耦合 法 能 最 准确 地 模拟 高 速 铅 芯 步枪 弹 侵 彻 钢 靶 
的 过 程 。FEM ЖІ EFG 法 受 限 于 网 格 ， 必 须 删 除 失效 单元 ， 因 此 仿真 过 程 中 质量 不 守恒 ， 
导致 侵 彻 阻力 偏 低 、 侵 彻 深度 偏 高 ， 且 不 能 准确 模拟 冲 塞 穿孔 过 程 。 

(2) 初始 光滑 长 度 参数 CSLH 为 1.4， 粒 子 近似 方法 为 重 整 化 近似 FORM 1， 失 效 
粒子 控制 方法 为 部 分 失效 、 不 保留 体积 响应 IEROD 1 的 参数 组 合 能 最 有 效 地 抑制 SPH 
方法 拉 伸 不 稳定 的 缺陷 ， 且 具有 较 好 的 计算 稳定 性 和 效率 。 

(3) 不 同 强度 钢 间 受 高 速 铅 心 弹 侵 彻 的 过 程 基 本 相同 ， 主 要 分 为 半 流 体 侵 彻 、 背 表 
面 损伤 、 剪 切 冲 塞 三 部 分 。 侵 彻 初 期 弹头 速度 较 高 ， 与 靶 板 相互 侵蚀 ; RUEIBZZAT 
其 抗 压强 度 ， 碍 体 受 力 面积 约 为 弹 体 头 部 截面 积 的 两 倍 。 侵 彻 过 程 中 ， 误 板 背 表面 受 整 
体 变形 影响 ， 主 要 承受 拉 应 力 ， 形 成 明显 隆起 并 发 生 拉 伸 损 伤 。 侵 彻 末 期 弹头 因 阻 力 而 
径 向 扩张 , 靶 板 受 侵蚀 而 厚度 减 小 , 形成 直径 与 弹 体 接 近 的 柱 形 剪 切 破坏 面 并 冲 出 塞 块 。 
中 低 强度 部 的 穿孔 以 半 流 体 侵 彻 为 主 ， 侵 彻 末 期 除 剪 切 冲 塞 以 外 伴 有 背 表面 拉 伸 破坏 ; 
高 强度 靶 板 的 穿孔 以 剪 切 冲 塞 为 主 ， 半 流体 侵 彻 识 度 相 对 较 低 且 与 抗 压 韧性 正 相 关 ， 剪 
切 破坏 面 较 明 显 ， 背 表面 仅 微 小 损伤 而 无 破坏 。 

(4) Rosenberg 半 流 体 侵 彻 一 Recht 剪 切 冲 塞 模型 可 定性 解释 铅 芯 弹 的 穿 甲 原理 ， 
但 受 限于 射 弹 结构 的 复杂 性 和 模型 自身 的 不 准确 性 ， 精 确 、 定 量 判断 击 穿 与 否 仍 依赖 于 
仿真 。 
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附录 A 文件 主要 内 容 


本 节 为 经 第 2 节 优 化 后 的 990 m/s M193 fU] 6 mm Armox 600 #8 k 文件 ,不 包含 过 
于 宛 长 的 节点 /粒子 CNODE) 、 单 元 (*ELEMENT SPH 和 *ELEMENT SOLID) 和 节 
点 集 (*SET NODE LIST) 关键 字 


*KEYWORD 
*CONTROL BULK VISCOSITY 
$# ql q2 type btype tstype 
1.5 0.06 -2 0 0 
*CONTROL CONTACT 
$  slsfac rwnpal islchk shlthk penopt thkchg orien enmass 
0.3 0.0 1 1 1 0 2 0 
$ usrstr usrfric nsbes interm xpene ssthk ecdt tiedprj 
0 0 0 0 0.0 0 0 0 
$ sfric dfric edc vfc th th sf pen sf unused 
0.0 0.0 0.0 0.0 0.0 0.0 0.0 
$ ignore frceng — skiprwg outseg | spotstp spotdel | spothin ^ unusedl 
2 1 0 1 0 1 0.5 
$ isym nserod rwgaps rwgdth rwksf icov swradf ithoff 
0 0 1 0.0 0.0 0 0.0 0 
$  shledg pstiff  ithcnt tdcnof ftall unused2 
1 0 0 0 1 0.0 
*CONTROL HOURGLASS 
$ ihq qh unused] 
5 0.05 
*CONTROL SPH 
$ ncbs boxid dt idim memory form start maxav 
1 01.00000E20 3 150 1 0.01.00000E18 
$ cont deriv ini ishow ierod icont javis isymp 
0 0 0 1 1 0 0 100 
$ ithk istab ql unused 


1 
*CONTROL_TERMINATION 


$  endtim endcyc dtmin endeng endmas unused 
5.00000E-5 0 0.001 0.0 100000.0 0 

*CONTROL TIMESTEP 

$  dtinit tssfac isdo tslimt dt2ms Ictm erode mslst 


0.0 0.6 0 0.0 0.0 0 1 0 
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$  dt2msf dt2mslc imscl unused 
0.0 0 0 0.0 
*DATABASE GLSTAT 
$# dt binary lcur ioopt 
6.00000E-8 3 0 1 
*DATABASE MATSUM 
$# dt binary lcur ioopt 
6.00000E-8 3 0 І 
*DATABASE BINARY D3PLOT 
$# dt lcdt beam npltc psetid 
5.00000Е-7 0 0 0 0 
$4 ioopt 
0 
*DATABASE EXTENT BINARY 
$ пеірһ neips maxint strflg sigflg epsflg ritflg engflg 
6 0 0 0 1 1 1 
$# cmpflg ieverp beamip dcomp shge stssz n3thdt ialemat 
0 0 0 1 1 1 2 
$# nintsld pkp sen sclp hydro msscl therm intout nodout 
0 0 0.0 0 2 0 
*BOUNDARY SPC SET 
$4 nsid cid dofx dofy dofz dofrx dofry dofrz 
1 0 1 1 1 1 1 
*CONTACT TIED NODES ТО SURFACE 
$# ssid msid sstyp mstyp sboxid mboxid spr mpr 
12 1 3 3 0 0 0 
$# fs fd dc vc vdc penchk bt dt 
0.0 0.0 0.0 0.0 10.0 0 0.01.00000E20 
$# sfs sfm sst mst sfst sfmt fsf vsf 
1.0 1.0 0.0 0.0 1.0 1.0 1.0 1.0 
$ Soft softscl Icidab maxpar sbopt depth bsort frcfrq 
2 0 0 0 3 5 0 
$ рептах tkhopt shlthk snlog isym i2d3d sldthk sldstf 
0 0 0 0 0 0 0 
*PART 
$# title 
Componentl'out 
$4 pid secid mid eosid hgid grav adpopt tmid 
1 1 1 1 0 0 0 
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*SECTION_SOLID 


aet 
0 
g e pr dtf 
79230.0 206000.0 0.3 0.0 
n с m tm 
0.175 0.00877 0.712 1537.85 
spall it dl d2 
2.0 0.0 -0.4 1.5 
erod efmin numint 
01.00000E-6 0.0 
sl s2 53 gamao 
1.49 0.0 0.0 1.93 
mid eosid hgid grav 
1 1 0 0 
hmin hmax sphini death 
0.2 3.0 0.01.00000E20 
mid eosid hgid grav 
3 3 0 0 
g е pr dtf 
43890.0 115000.0 0.31 0.0 
n с m tm 
0.42 0.009 1.68 915.85 
spall It 41 а2 
2.0 0.0 0.0 2.65 
erod efmin numint 
01.00000E-6 0.0 


$4 весій elform 
1 1 
*MAT JOHNSON COOK 
$4 mid ro 
17.85000E-9 
$4 a b 
1580.0 958.0 
$4 cp pc 
4.520000E8 0.0 
$4 d5 c2/p 
0.0 0.0 
*EOS GRUNEISEN 
$4  eosid с 
] 4570000 
$# v0 
0.0 
*PART 
$# 
in] 
$# pid secid 
12 5 
*SECTION SPH 
$ весі cslh 
5 1.4 
*PART 
$# 
Jacket 
$# pid secid 
13 5 
«МАТ JOHNSON СООК 
$# mid ro 
38.80000E-9 
$# a b 
111.69 504.69 
$# cp pe 
3.760000E8 0.0 
$# d5 c2/p 
0.0 0.0 
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vp rateop 
0.0 0.0 
tr epso 
22.0 1.0 
d3 d4 
-0.5 0.0 
a е0 
0.0 0.0 
title 
adpopt аша 
0 
start sphkern 
0.0 1 
title 
adpopt аша 
0 
vp rateop 
0.0 0.0 
tr epso 
22.05.00000E-4 
d3 d4 
-0.62 0.28 


*EOS_GRUNEISEN 
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$4 сова с sl s2 s3 gamao 
3 3940000 1.49 0.0 0.0 1.97 
$# v0 
0.0 
*PART 
$# 
core 
$# pid secid mid eosid hgid grav 
14 5 4 4 0 0 
«МАТ JOHNSON COOK 
$# mid ro g e pr dtf 
41.13400E-8 7000.0 20090.0 0.435 0.0 
$# a b n с т tm 
1.0 63.0 0.163 0.41 1 326.85 
5 ср рс spall it dl d2 
1.690000Е8 0.0 2.0 0.0 2.0 1.77 
$ d5 c2p 
0.0 0.0 01.00000E-6 0.0 
*EOS GRUNEISEN 
$4 сова C sl s2 s3 gamao 
4 2006000 1.429 0.0 0.0 2.74 
$# v0 
0.0 
*INITIAL VELOCITY GENERATION 
$#nsid/pid styp omega VX Vy VZ 
13 2 19314.8 0.0 0.0 990000.0 
$# xc yc zc nx ny nz 
0.0 0.0 0.0 0.0 0.0 0.0 
$ sid styp omega VX Vy VZ 
14 2 19314.8 0.0 0.0 990000.0 
5 хс ус ZC nx ny nz 
0.0 0.0 0.0 0.0 0.0 0.0 
*DEFINE COORDINATE SYSTEM 
$# cid xo yo zo хі уі 
93 0.0 0.0 0.0 1.0 0.0 
$# xp yp zp 
0.0 1.0 0.0 


a e0 
0.0 0.0 
title 
adpopt tmid 
0 
vp rateop 
0.0 0.0 
tr epso 
22.05.00000E-4 
d3 d4 
-3.4 0.0 
unused 
a e0 
0.0 0.0 
ivatn icid 
0 93 
phase irigid 
0 0 
ivatn icid 
0 93 
phase unused2 
0 0 
zl cidl 
0.0 0 
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修改 记录 


一 、 毕 业 论文 〈 设 计 ) 内 容重 要 修改 记录 
指导 教师 : 

1) 扩充 理论 建 模 章 

2) 调整 附录 A 格式 
评阅 教师 : 

1) 修改 公式 缩 进 格式 

2) 调整 图 表格 式 
答辩 委员 会 : 
1) 扩充 摘要 
2) 修改 页 眉 格式 
二 、 毕 业 论文 〈 设 计 ) 外 文 翻译 修改 记录 
1) 修改 图 表 标 题 格式 

2) 调整 标题 分 级 格式 
三 、 毕 业 论文 〈 设 计 ) 正式 检测 重复 比 
总 文字 复制 比 1.7%， 重 复 字 数 529 


记录 人 (签字 ) . FRM 
指导 教师 сөз. A A £, 
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Spo й 


本 论文 是 在 大 连理 工大 学 机 械 工程 学 院 薄 其 乐 工程 师 和 工程 力学 系 韩 啸 副教授 的 
指导 下 悉心 完成 的 。 在 论文 完成 之 际 ， 这 里 有 一 份 简短 的 致谢 名 单 : 

首先 ， 要 感谢 大 连理 工大 学 机 械 工程 学 院 四 年 来 对 我 的 培养 。 没 有 学 校 的 支持， 学 
业 论 文 的 完成 是 难以 想象 的 。 

其 次 ， 最 要 感谢 的 是 工程 力学 系 韩 啸 老师 在 研究 工作 上 的 悉心 指导 。 在 这 里 ， 我 获 
得 了 作为 一 名 普通 的 本 科 生 本 不 可 能 取得 的 学 术 成 就 。 不 仅 收 获 了 科研 上 的 技术 知识 ， 
也 学 到 了 做 研究 的 思路 和 方法 。 

随后 ， 要 感谢 两 年 半 来 工程 力学 系 李 晓 杰 教授 对 我 的 关怀 和 帮助 ， 在 初 入 本 研究 领 
域 时 为 我 提供 基础 的 技术 指导 和 研究 方向 的 指引 。 

最 后 ， 感 谢 北 京 市 系 源 学 校 秦 去 和 赵 晶 老师 的 长 期 支持 与 鼓励 。 

愿 未 来 能 为 高 速 侵 彻 领域 的 研究 贡献 自己 的 一 份 力量 。 


